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RESUMO 


A necessidade de novas tecnologias em nossa sociedade fez com que houvesse um maior 
interesse em pesquisas de materiais utilizados em aplicações em eletrônica, isso tornou os 
materiais cristalinos multiferróicos muito promissores e de grande interesse acadêmico e 
tecnológico devido a correlação entre as propriedades ferroelétricas e (anti)ferromagnéticas. 
Podemos destacar entre estes materiais os compostos com estruturas do tipo perovskita, 
bem como os que surgem a partir de dopagens no BiFeO,. No presente trabalho foram 
processados os pós cerâmicos do composto (Bi, Nd,)(Fe,,Co JO, com (x=0,0; 0,1 e 0,2) e 
(y=0,0; 0,01 e 0,02), obtidas por moagem em altas energias, seguido de sinterização em 
atmosfera livre. Também foi utilizada a técnica de sinterização rápida seguida de choque 
térmico (resfriamento rápido). Fazendo uso da técnica de difração de raios X, e do método 
de refinamento estrutural de Rietveld concluiu-se que as composições possuem simetria 
romboédrica com grupo espacial R3c. Os cálculos de densidade eletrônica permitiram uma 
melhor análise quanto ao caráter de ligação entre os íons de Biº*/Ndº%* e OZ, e entre o 
íons de Fe*t/Co** e O? do composto (Bi-Nd)(Fe,-,Co)O.. Pudemos observar que com 
a adição de cobalto as ligações entre os íons de ferro e oxigênio tendem a se tornar mais 
covalentes. O mesmo não aconteceu com as ligações entre os íons de bismuto e oxigênio. 
Para as amostras dopadas com neodímio, pudemos observar que com x = 0,1 houve a 
total supressão da ligação covalente (Lone Pairs) entre os íons de bismuto e oxigênio e um 
aumento significativo na covalência entre os íons de ferro e oxigênio. 


PALAVRAS-CHAVE: Multiferróicos, Refinamento Rietveld, Estrutura Eletrônica, MEM. 
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ABSTRACT 


The need for new technologies in our society caused a greater interest in research on 
materials used in electronics applications, it made the multiferroic crystalline materials very 
promising and of great and technological academic interest materials due to correlation 
between the ferroelectric and (anti)ferromagnetic properties. Among these materials, 
compounds with perovskite structures, as those arising from the doping of BiFeO, can be 
mentioned. In this work ceramic powders of (Bi, Nd,)(Fe, Co JO, compositions, for (x = 0.0, 
0.1 and 0.2) and (y = 0.0, 0.01 and 0.02), were processed using high-energy ball- milling 
followed by sintering in a free atmosphere. Also it was used the fast sintering technique 
followed by quenching. By using the X-ray diffraction technique, and structural refinement 
via Rietveld method it is concluded that all compositions have rhombohedral symmetry 
with space group R3c. The electron density calculations allowed better analysis on the 
bonding character between the Biºt/Ndº* and O? ions, and between the Fe**/Co** and O? 
ions in the (Bi, «Nd (Fe, Co )O3 compositions. We could observe that with the addition 
of cobalt ions the bonding between iron and oxygen ions tend to become more covalent. 
This did not happen with the bonding between bismuth and oxygen ions. For neodymium 
doped samples, we observed for x = 0.1 the complete disappearance of the covalent bond 
(Lone Pairs) between the bismuth and oxygen ions and a significant increase in covalency 
between iron and oxygen ions. 


KEYWORDS: Multiferroic, Rietveld Refinement, Electronic Structure, MEM. 
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INTRODUÇÃO 


A necessidade de novas tecnologias em nossa sociedade fez com que houvesse um 
maior interesse em pesquisas de materiais utilizados em aplicações com eletrônica, como 
sensores e atuadores. Conhecidos como ferróicos primários os materiais com propriedades 
ferroelétricas, ferromagnéticas e ferroelásticas, destacam-se em pesquisas e aplicações 
tecnológicas. Em geral esses materiais respondem a somente um tipo de excitação externa 
(campos elétricos, campos magnéticos ou tensões mecânicas) produzindo uma resposta 
única [1], [2]. 

As tentativas de combinar em um único sistema, as propriedades ferromagnética 
e ferroelétrica começaram em 1960 [3], na União Soviética. Devido à combinação dessas 
propriedades ferróicas, mais tarde estes materiais foram chamados de “multiferróicos”. Por 
algum tempo, esse campo de pesquisa ficou inerte sendo ainda desconhecido. O interesse 
por estes materiais teve um aumento por volta de 2001-2003, devido a diversos fatores, 
especialmente, à melhoria da técnica de preparação e estudo de filmes finos de óxidos. 
Tal avanço permitiu a produção de bons filmes finos de materiais ferroelétricos, e abriu 
a possibilidade de utilizar estes sistemas, por exemplo para a produção de memórias 
ferroelétricas; e a descoberta de vários novos sistemas multiferróicos com propriedades 
interessantes, em especial um forte acoplamento entre graus de liberdade ferroelétricos 
e magnéticos [3]. Atualmente é possível a produção de novos materiais multiferróicos 
multifuncionais nos quais coexistem duas ou mais propriedades ferróicas primárias, 
por exemplo, a coexistência das propriedades magnética e ferroelétrica, chamada de 
magnetoelétrica. O efeito magnetoelétrico pode ser definido como um acoplamento 
de respostas aos estímulos externos, na qual um campo elétrico pode produzir uma 
magnetização ou um campo magnético pode produzir uma polarização elétrica [1], [2]. 

Nas últimas décadas a investigação sobre as propriedades físicas destes materiais 
mostraram que a multifuncionalidade, relacionada com a coexistência e acoplamento entre 
as propriedades ferróicas primárias, é muito pequena ou basicamente proibida por razões 
de simetria. Por exemplo, nos materiais magnéticos e ferroelétricos, as propriedades 
magnéticas são provenientes das camadas d incompletas enquanto que a ferroeletricidade 
surge das camadas d vazias, respectivamente, isso torna a coexistência de (anti) 
ferromagnetismo e ferroeletricidade rara [2]. Tanto no ordenamento magnético como no 
ferroelétrico a forma como as camadas eletrônicas estão preenchidas é fundamental para 
que ocorra um estado ordenado. 

O aumento de novas classes de materiais heterogêneos e compósitos nos últimos 
anos permitiu novos estudos frente as propriedades magnéticas e ferroelétricas [3]. Podemos 
destacar alguns materiais, como o BiMnO, e o BiFeO,, que apresentam íons magnéticos 
(Mnº+ e Fe%+) e também são ferroelétricos [4]. Nesses materiais a ferroeletricidade se deve 
ao íon de Bi, que possui dois elétrons no orbital 6s, cnamados de “Lone Pair”, que se 


Introdução 


desloca do centro de simetria [4]. O estudo das ligações químicas e sua relação com as 
propriedades ferróicas é de grande importância para se entender a origem das propriedades 
multifuncionais apresentadas por materiais baseados no BiFeO,. 

Portanto, devido à importância dos multiferróicos na compreensão da ciência e 
em aplicações tecnológicas, o presente trabalho realizou estudos acerca de compostos 
baseados no BiFeO, com a substituição de óxidos de neodímio e cobalto. 


11 OBJETIVOS 
Os objetivos principais desta dissertação são: 


* As caracterizações estruturais das amostras do composto (Bi, Nd (Fe, ,Co,) 
O.. 
*- O estudo da ligação química e sua possível relação com as propriedades fer- 

róicas. 

Para o estudo das estruturas cristalinas do composto (Bi Nd (Fe, yCo)O, com (x 
=0,0;0,1€e0,2)e (y = 0,0; 0,01 e 0,02), fizemos uso da difração de raios X. Em seguida 
os parâmetros estruturais foram refinados por método de Rietveld utilizando os programas 
Rietan' e FullProf?. 

Ainda, realizamos um estudo sobre as implicações de nossos resultados no que se 
refere às ligações químicas e suas forças. Os estudos utilizando as densidades eletrônicas 
do composto (Bi, Nd)(Fe,,Co)O, foram realizados por meio do Método da Máxima 
Entropia calculado através do programa Dysnomiaº. Os mapas das densidades eletrônicas 
foram visualizados no programa VESTA.. O cálculo das intensidades medidas pela mínima 
densidade eletrônica (MED) tornou possível as análises das ligações químicas. 


21 ORGANIZAÇÃO DA DISSERTAÇÃO 


Inicialmente no capítulo 2, apresentamos a fundamentação teórica de nosso 
trabalho, fazendo uma revisão da literatura a respeito da ferroeletricidade, magnetismo, 
materiais multiferróicos e ainda uma revisão bibliográfica a respeitodo BiFeO, e dopagens 
realizadas com este composto. 

Em seguida, no capítulo 3 são apresentadas técnicas de preparação de amostras, 
como moagem e sinterizações, e ainda a técnica de caracterização estrutural por difração 
de raios X. Também são apresentados os métodos de refinamento estrutural, cálculos de 
densidade eletrônica e análises quanto as ligações químicas. 

No capítulo 4 são descritos os procedimentos experimentais utilizados neste 
trabalho, o método de preparação das amostras, as caracterizações, os refinamentos e os 


1 Programa gratuito, disponível em: http://fujioizumi.verse.jp/download/download Eng.html. 
2 Programa gratuito, disponível em: https://www.ill.eu/sites/fullprof/php/downloads.html. 

3 Programa gratuito, disponível em: http://jp-minerals.org/dysnomia/en/. 

4 Programa gratuito, disponível em: http://jp-minerals.org/vesta/en/download.html. 


Introdução 


cálculos da densidade eletrônica. 
E por fim a discussão dos resultados obtidos nesta pesquisa são apresentados no 
capítulo 5. 
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FUNDAMENTAÇÃO TEÓRICA 


11 MATERIAIS MAGNÉTICOS 


Materiais magnéticos são aqueles cujos átomos apresentam um momento magnético. 
Esses materiais apresentam uma transição de fase de uma temperatura mais elevada, na 
qual esses momentos não possuem um alinhamento, para uma fase a temperatura mais 
baixa na qual ocorre o alinhamento desses momentos magnéticos. Esse alinhamento é 
criado por interações de troca quânticas, para as quais a energia magnética associada pode 
favorecer um alinhamento paralelo, antiparalelo ou algumas variações de um alinhamento 
antiparalelo [5]. Nos materiais ferromagnéticos, que exibem magnetização espontânea à 
temperatura ambiente e mesmo na ausência de campo magnético externo possui momento 
magnético (mM 0) , encontramos um alinhamento paralelo e de mesma magnitude dos 
momentos magnéticos dentro dos domínios (Figura 2-1b). Acima da temperatura crítica T.. 
(temperatura de Curie), a magnetização espontânea desaparece e os spins se rearranjam 
dando origem a um material paramagnético (Figura 2-1a). 

Na Figura 2-1c, temos uma orientação ferrimagnética, nesta o material possui 
momentos de dipolo permanentes que interagem entre si, causando alinhamento 
antiparalelo e desigual das sub-redes magnéticas (m 0) . 

Na Figura 2-1d, observamos um ordenamento antiferromagnético, no qual os spins 
são dispostos de forma antiparalela com momento total nulo (Mm = 0) em temperaturas 
inferiores a T,. A ordem antiferromagnética desaparece acima da temperatura crítica (T,), 
e o material torna-se paramagnético [6]. 


ei um 
bu. 
“00 dóddo dIdIA dedos 


Paramagnético Ferromagnético Ferrimagnético Antiferromagnético 


Figura 2-1:a) Orientação paramagnética, b) orientação ferromagnética, c) orientação ferrimagnética, d) 
orientação antiferromagnética. 


Ao analisarmos uma curva de histerese conhecemos informações referentes às 


propriedades magnéticas de um material. A Figura 2-2 ilustra um ciclo de histerese que 
relaciona a magnetização (M) e o campo magnético (H). 
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Figura 2-2: Ciclo típico de histerese M-H para materiais ferromagnéticos [1]. 


O ciclo é gerado a partir da magnetização medida em função do campo magnético 
aplicado. No ponto “M ” a magnetização atinge a saturação, praticamente todos os domínios 
magnéticos estão alinhados. Quando o campo magnético H é reduzido a zero, a curva 
move-se a partir do ponto de saturação “M ” para o ponto de magnetização remanescente 
“M/. O campo magnético necessário para remover a magnetização residual do material é 
chamado de campo magnético coercivo (H ). Se aumentarmos o campo coercitivo, teremos 


uma magnetização com orientação contrária [7]. 


21 MATERIAIS FERROELÉTRICOS 

Os materiais ferroelétricos diferenciam-se dos demais dielétricos por apresentar 
uma polarização elétrica espontânea que pode ser revertida com a aplicação de um campo 
elétrico externo. Dessa forma exibem um efeito de histerese entre a polarização e o campo 
elétrico externo aplicado. 

No estado ferroelétrico, os centros das cargas positivas não coincidem com o centro 
das cargas negativas, o desequilíbrio de cargas dá origem a momentos de dipolo, em casos 
especiais, é possível inverter a direção dos momentos de dipolo, neste caso temos um 
material ferroelétrico. Ou seja, os materiais ferroelétricos são “dielétricos especiais”, por 
possuírem uma polarização espontânea mesmo na ausência de campo elétrico externo. A 
orientação da polarização espontânea pode ser modificada pela aplicação de um campo 
elétrico externo. Para entendermos melhor consideremos a cela unitária de estrutura 
perovskita (vide Figura 2-3), típica do titanato de bário (BaTiO,), na fase ferroelétrica (288 
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K<T<400K) [6], [7]. 


OB CTDO 


Figura 2-3: Estrutura perovskita típica do titanato de bário (BaTiO,) ; e segundo a vista lateral de uma 
das faces, a qual ilustra o deslocamento do centro da face [10]. 


Observa-se que os íons Ba” estão localizados nos vértices da célula unitária, que 
apresenta simetria tetragonal. O momento de dipolo surge dos deslocamentos relativos 
dos íons o” e Tm das suas posições simétricas. Os íons de oxigênios estão localizados, 
ligeiramente abaixo, dos centros de cada uma das seis faces, enquanto o íon Ti“ encontra- 
se localizado logo acima do centro da célula unitária. Com isso temos um momento de 
dipolo iônico permanente associado à cada célula unitária. Mas se a temperatura estiver 
acima da temperatura de Curie (-400 K) a célula unitária se torna cúbica e os fons assumem 
posições simétricas fazendo com que o comportamento ferroelétrico deixe de existir (Figura 
2-4). 


(a) 


Figura 2-4: A estrutura do cristal de titanato de bário (estrutura de tipo perovskita). (a) Acima da 
temperatura de Curie a célula é cúbica; (b) abaixo da temperatura de Curie a estrutura é tetragonal com 
íons de Ba eTi deslocados em relação aos íons de O [10]. 
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O comportamento ferroelétrico é observado em intervalos de temperatura 
delimitados pela Temperatura de Curie (T,), acima desta temperatura a ordem ferroelétrica 
é destruída e o material assume uma fase não polarizada, denominada fase paraelétrica. A 
propriedade ferroelétrica de um material surge devido a uma pequena distorção da estrutura 
cristalográfica de sua fase paraelétrica, de maneira que a simetria na fase ferroelétrica é 
menor que na fase paraelétrica. Sendo assim, pode-se dizer que o ponto de Curie é a 
temperatura a partir da qual a estrutura cristalina do material se torna mais simétrica [7], 
[9]. A resposta dos domínios ferroelétricos frente a um campo elétrico externo aplicado 
em um ferroelétrico é representada pela curva denominada ciclo de histerese ferroelétrica 
(Figura 2-5). 

Assim como os materiais ferromagnéticos, os ferroelétricos possuem um ciclo de 
histerese (vide Figura 2-5), caracterizado pelos pontos de polarização de saturação (P), 
polarização remanescente (P)) e campo elétrico coercitivo (E ). A polarização remanescente 
se deve à orientação dos domínios, que permanece no material na ausência de um 
campo elétrico externo, sendo invertida somente sob a aplicação de um campo elétrico de 
amplitude superior à amplitude do campo coercitivo. 


Figura 2-5: Ciclo típico de histerese P-E em ferroelétricos [10]. 


Os domínios podem experimentar uma nova reorientação se o campo é revertido 


pelo segmento GHC, e o ciclo de histerese se repete [9]. 


31 MATERIAIS MULTIFERRÓICOS 


São denominados multiferróicos os materiais que apresentam duas, ou mais das 
chamadas ordens ferróicas, ou seja, ferroeletricidade, ferromagnetismo e ferroelasticidade 


Fundamentação Teórica 


[4], [12], [13]. 

Nos materiais multiferróicos, tem-se, por exemplo, a possibilidade de controlar a 
magnetização de um material com um campo elétrico ou a polarização elétrica com um 
campo magnético (efeito magnetoelétrico) [14]. Se o composto apresenta correlação entre 
as propriedades mecânicas e elétricas é chamado de piezoelétrico, pois uma tensão no 
cristal pode gerar um campo elétrico. Se o composto possui correlação entre propriedades 
mecânicas e magnéticas, este é conhecido como piezomagnético, pois uma mudança nas 
dimensões do cristal pode gerar uma magnetização espontânea [13]. 

Podemos ver a representação destas interações na Figura 2-6, na qual P representa 
a polarização elétrica; E representa o campo elétrico aplicado; M representa a magnetização 
do material; H representa o campo magnético aplicado; £ representa o estresse mecânico e 


o representa a tensão mecânica aplicada [15]. 


Figura 2-6: Diagrama esquemático do controle de fases em materiais ferróicos e multiferróicos. Nos 
materiais ferróicos, o campo elétrico E controla a polarização P (em amarelo), o campo magnético H 
controla a magnetização M (em azul) e uma pressão ou tensão mecânica o (comumente chamada de 
“stress”) controla o “strain” (deformação relativa) e, que é compactação ou contração (em vermelho) 
[16]. 


Entre as interações vistas na figura 2-6, a magnetoelétrica tem atraído grande 
interesse em pesquisas por apresentar o acoplamento entre os parâmetros de ordem 
magnética e elétrica. A dificuldade de se encontrar materiais com esse tipo de acoplamento 
se deve ao número reduzido de compostos com estas características. Estudos teóricos 
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mostram que, em escala atômica, os mecanismos que conduzem às ordens magnéticas e 
ferroelétrica são exclusivos, pois requerem orbitais vazios ou parcialmente preenchidos das 
camadas atômicas dos metais de transição [17]. Dentre os raros materiais que possuem 
propriedades multiferróicas, destacamos algumas perovskitas do tipo ABO,, tais como 
HoMnoO, e BiFeO,. 


41 ESTRUTURA PEROVSKITA 


Alguns cristais ferroelétricos importantes podem apresentar uma estrutura conhecida 
como perovskita, genericamente esta classe de estrutura apresenta a fórmula química 
ABO,, onde o íon A tem 12 átomos de coordenação e fica no vértice de um cubo ideal; o 
íon B (menor que o íon A) tem 6 átomos de coordenação e fica no centro do cubo; os íons 
de oxigênio ficam na face do cubo como mostrado a Figura 2-7a. Na Figura 2-7b a mesma 
estrutura pode ser visualizada sob o ponto de vista dos octaedros de oxigênio BO,. Nesta 
configuração a célula unitária é representada com os cátions menores nos vértices de um 
cubo, os cátions maiores no centro da estrutura e os íons de oxigênio nas extremidades 
[8], [9], [18]. 


Figura 2-7: (a) típica estrutura perovskita ABO, e (b) estrutura perovskita do ponto de vista dos 
octaedros de oxigênio BO, [18]. 


Para que a estrutura perovskita permaneça estável é necessário que exista um 
arranjo BO, estável e que o cátion A tenha um tamanho adequado para ocupar o interstício 
gerado pelos octaedros. Neste sentido, Goldschmidt [19], efetuou estudos ainda na década 
de 1920 e estabeleceu um fator de tolerância que é dado pela expressão: 


(r,+ro) 


V2(15 +19) na 


f= 
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Onde ra, rB e ro são os raios iônicos dos elementos que ocupam os sítios A, Be O, 
respectivamente, na estrutura perovskita. Um átomo é considerado estável numa 
estrutura perovskita cúbica quando este fator de tolerância está próximo de um. No 
entanto, a estrutura cúbica é apenas observada na fase paraelétrica dos materiais 
ferroelétricos de estrutura perovskita [19]. 

Em geral, a maior parte das perovskitas ABO, apresenta perda de simetria para 
baixas temperaturas devido a distorções provocadas por deslocamentos atômicos, 
nas quais o íon central é deslocado da posição central e a célula unitária se deforma 
assumindo uma célula em que não há um centro simétrico, como tetragonal, romboédrica 


ou monoclínica, vide Figura 2-8. 


Romboédrica Monoclínica Tetragonal 
a=b=c atbéc az=béc az=b=c 
a=B=Y *90º B=Y=90ºa ao=B=Y az=p=Y=90º 


Figura 2-8: Esboço das transições estruturais que ocorrem com a estrutura perovskita com o aumento 
da temperatura da esquerda para direita, adaptado de [5]. 


Contudo, essa perda de simetria provoca um rearranjo que irá proporcionar o 
fenômeno da ferroeletricidade. 

Embora exista um grande número de perovskitas simples ABO,, o número de 
compostos aumenta significativamente quando um ou mais dos íons originais são 
substituídos por outros íons. Na maioria dos casos essa substituição acontece nos sítios 
dos cátions e gera um grupo enorme de compostos conhecido como perovskitas compostas, 
com fórmula química AN'BB'O, [9]. 


51 REVISÃO BIBLIOGRÁFICA 


5.1 BiFeO3 


O BiFe0O, é provavelmente o material multiferróico mais promissor, e certamente 
o mais estudado que possui como base uma perovskita ferroelétrica [20]. O estudo do 
BiFeO,, como um material multiferróico importante, foi iniciado em 1960 por Smolenskii 
e colaboradores [21]. Desde então, o número de trabalhos publicados apresentou um 
aumento discreto nas publicações até 2002, mas a partir desse ano o número de trabalhos 
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sobre esse material teve um aumento significativo, devido ao seu alto potencial tecnológico. 
Espera-se o uso deste material em spintrônica [22, 23], no desenvolvimento de memórias 
ferroelétricas [24] e de capacitores multicamadas [25]. 

A ferrita de bismuto (BiFeO,) é um dos materiais multiferróicos, no qual parâmetros 
estruturais de ordens elétricas e magnéticas são acoplados um ao outro, resultando 
em propriedades ferróicas simultâneas, como ferroeletricidade, antiferromagnetismo e 
ferroelasticidade. Diferentemente da maioria dos outros multiferróicos, BiFeO, apresenta 
propriedades multiferróicas à temperatura ambiente devido a sua ferroeletricidade, possui 
temperatura de Curie (T, = 820ºC) e antiferromagnetismo com uma temperatura de Néel (T, 
= 380ºC). Esses fatores tornam o BiFeO, um material de grande potencial para aplicações 
industriais [26]. 

No entanto, este material apresenta vários problemas nas propriedades físicas 
primárias. Sua aplicação ainda não é viável devido às dificuldades em restringir a formação 
de fases de impurezas e de poros quando desenvolvidos em bulk, assim como aliviar as 
altas densidades de corrente de fuga, que levam à falha dielétrica do material [27]. As fases 
secundárias que ocorrem mais frequentemente durante as sínteses do BiFeO, são as fases 
Bi,Fe,0, e a Bi, FeO,, [3]. 


5.2 Estrutura e Polarização do BiFeO, 


O BiFeO, possui uma estrutura perovskita romboédrica distorcida (grupo espacial 
R3c), na sua fase ferroelétrica pode ser representado por uma estrutura hexagonal, como 
podemos ver na Figura 2-9. 


Figura 2-9: Célula unitária hexagonal do BiFeO,, construída com o programa Vesta a partir dos dados 
da referência [28]. 
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Os parâmetros de rede e as posições atômicas do BiFeO, romboédrico (R3c), retirados 


da referência [28], estão apresentados na Tabela 2-1. 


Grupo Espacial R3c 
a=b=5,587 À c=13,867 À 
íon x y z 
BiR+ 0 0 0 
Fe 0 0 0,2212 
(op 0,443 0,012 0,9543 


Tabela 2-1: Dados estruturais do BiFeO, para o grupo espacial R3c, adaptado da referência [28]. 


A Figura 2-10 mostra uma célula unitária do tipo Perovskita do composto BiFeO, 


construída com os dados obtidos da referência [28]. 


b 


Figura 2-10: Célula unitária tipo Perovskita do composto BiFeO,. Construída com o programa VESTA a 
partir dos dados da referência [28]. 


No programa VESTA, foi possível calcular as distâncias de ligação entre os íons nos 
planos (001), (200) e (002), estas distâncias são apresentadas na Tabela 2-2. 
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Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (A) Ligação (À) 

06 — Bit d7 3,400 Fe — O2 d6 1,940 

06 — Bi2 d2 2,523 Fe-— 03 dz 1,940 

06 — Bi3 [oh] 2,310 Fe-—- 04 dg 2,120 

06 — Bi4 d4 3,210 Fe-— 05 dg 2,120 

Fe-01 ds 1,940 Fe — 06 dio 2,120 


Tabela 2-2: Distâncias de ligação obtidas partir dos dados da referência [28] e calculadas com o 
programa VESTA. 


A estrutura R3c não é centrossimétrica, isso permite deslocamentos relativos ao 
longo do plano (111), de maneira especial o Bi em relação ao octaedro FeO, (Figura 2-11). 


Figura 2-11: Esquema de uma estrutura R3c romboédrica, onde os cátions são deslocados ao longo da 
direção [111] relativo aos ânions, e os octaedros de oxigênio giram com sentido contrário em torno da 
linha central de [111]. 


Cálculos da polarização espontânea do BiFeO,, com estrutura R3c romboédrica, 
ao longo da direção [111] sugerem um valor entre 90 e 100 uC/cm', devido aos Lone Pairs 
(pares de elétrons não compartilhados) que influenciam nas interações de ligação que 
ocorrem entre os sub-níveis dos orbitais 6s” dos íons de Bi” e 2p” dos íons de O” que 
provocam uma distorção na geometria do octaedro de FeO, [29]. 


5.3 Efeitos estruturais com a dopagem do composto BiFe0O, 


Um dos problemas do BiFeO, referido por diversos autores é a baixa resistividade 

da mostra à temperatura ambiente, tornando difícil a observação do ciclo de histerese. 
3+ 3+ 3+ 3+ 

A substituição parcial de Bi por íons terras-raras (TR) como (La , Nd , Gd , Dy ) entre 
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outros, é conhecida por induzir uma resposta ferromagnética, diminuindo a condução 
e eliminando as lacunas de oxigênio que atuam como portadores de carga que geram 
condução do cristal [30], [31]. 

Mathe e colaboradores sintetizaram amostras policristalinas de Bi, ,Nd,FeO, através 
do método de reação de estado sólido. Nos estudos de difração de raios X, verificaram 
que a estrutura cristalina da cerâmica para x < 0,2 foi descrita pela simetria romboédrica 
enquanto que para x > 0,3 por simetria ortorrômbica, no entanto para x = 0,2, o material 
cristalizava numa estrutura triclínica [32]. 

Zhang e colaboradores analisaram o Bi, R,FeO, (xvariando de0a 1,R=La, Nd, Sm, 
Eu e Tb), utilizaram a sinterização de estado sólido convencional em pressão atmosférica. 
As difrações de raios X foram refinadas pelo método de Rietveld, onde se verificou que 
a estrutura cristalina do BiFeO, se transforma de romboédrica para ortorrômbica com o 
aumento do teor das diferentes terras raras. Observaram que a estrutura do Bi, Nd,FeoO, 
é romboédrica (R3c) para x < 0,2 e que apresenta uma transição estrutural para a simetria 
ortorrômbica (Pnma) para 0,2 <x < 0,4, acima de 0,4 até x = 1 a estrutura continua a 
ser ortorrômbica (Pnma). Pode se verificar também que há uma diminuição no volume da 
célula unitária com o aumento da dopagem [33]. 

Kumar e Varshney analisaram amostras policristalinas de Bi, ,Nd,FeoO, (x=0,0; 0,175; 
0,20) preparadas pelo método de co-precipitação química. O efeito da dopagem sobre as 
propriedades estruturais do material foi estudado utilizando a difração de raios X (DRX), 
fazendo o refinamento estrutural pelo método de Rietveld, verificaram que para BiFeO, a 
estrutura é romboédrica (R3c) e após a dopagem com Nd (0,175 e 0,20) a estrutura do 
material é triclínica (P1) [23]. 

Kumar e colaboradores prepararam amostras policristalinas de Big RçpFeo, 
(TR = La, Nd, Dy) e refinaram os difratogramas de DRX utilizando o método de Rietveld, 
verificaram que a amostra BLFO (BiFeO, dopado com La +) possui estrutura romboédrica 
(R3c), já as amostras BNFO e BDFO revelaram as seguintes mudanças de fase R3c — P1 
e R3c — Pnma respectivamente [34]. 

Dzik e colaboradores estudaram o Bi, Nd,FeO, para x = 0,1-0,4. O método de óxido 
misto seguido por sinterização sem pressão foi utilizado para a fabricação das cerâmicas. 
Verificaram que a concentração de neodímio tem um efeito sobre o tamanho das partículas, 
um aumento no teor de neodímio causou uma diminuição no tamanho médio dos grãos da 
cerâmica. Verificou-se que a amostra de Bi, Nd,FeO, com x s 0,2 deve ser descrita por 
simetria romboédrica, grupo espacial R3m (nº 160), e para x > 0,3 a estrutura cristalina 
do Bi, Nd,FeO, deve ser descrita por simetria ortorrômbica, grupo espacial Pnma (nº 62), 
havendo uma transição de fases para 0,2 <x< 0,3 [35]. 

Em estudos de filmes finos preparados com Bi, Nd,FeO, (x = 0,06-0,20), Xue 
e colaboradores observaram que para x < 0,20 a estrutura cristalina possui simetria 
romboédrica, mas para x = 0,20 foi observado uma nova simetria estrutural, sendo esta 
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ortorrômbica [36]. 
Na Tabela 2-1 apresenta de uma forma resumida as diferentes propostas de estrutura 
para o composto Bi Nd,FeO.. 


Dopagem Estrutura proposta Referências 

x<0,2 Romboédrica Mathe et al [32] 
x=0,2 Triclínica 

x20,3 Ortorrômbica 

x<0,2 Romboédrica, R3c Zhang et al [33] 
0,2=x<0,4 R3c—>Pnma 

0,4<x=<1 Ortorrômbica, Pnma 

x=0,0 Romboédrica, R3c Kumar e Varshney [23] 
x=0,175 e 0,20 Triclínica, P1 

x=0,0 Romboédrica, R3c Kumar, Sharma e Varshney [34] 
x=0,20 Triclínica, P1 

x<0,2 Romboédrica, R3m Dzik et al [35] 
0,2<x<0,3 R3c—>Pnma 

Xz0,3 Ortorrômbica, Pnma 

x<0,20 Romboédrica Xue et al [36] 

x=0,2 Ortorrômbica 


Tabela 2-3: Diferentes propostas de estruturas para Bi, ,Nd,Feo, 


Luo e colaboradores analisaram a amostra de Bi(FO, 260, 0/04 utilizando difração 
de raios X, e encontraram que os picos de difração são semelhantes aos do BiFeO, puro, 
indicando que nenhuma transformação estrutural foi induzida, ou seja, a dopagem com 
Cobalto apresenta estrutura romboédrica R3c [37]. 

Em análise aos efeitos da dopagem no composto (Bi, Nd,)(Fe,s,/C0,,;)03 com (x = 
0,06; 0,09; 0,12 e 0,15), Xue e colaboradores, utilizando difração de raios X e refinamento 
estrutural de Rietveld, mostraram que houve uma transição de uma simetria romboédrica 


para triclínica com o aumento na dopagem com neodímio [38]. 
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TÉCNICAS E MÉTODOS 


11 MOAGEM EM ALTAS ENERGIAS 


A moagem em altas energias é uma técnica desenvolvida no final dos anos 
sessenta pela “International Nickel Company” [39], esta técnica de processamento permite 
a produção de materiais homogêneos a partir da mistura de pós-precursores, reunidos 
com esferas de aço ou outro material de alta dureza, inseridos em um vaso de moagem, 
em movimento energético. Com vibração ou rotação as esferas colidem com as paredes 
do vaso resultando em uma prensagem do pó a cada impacto. No decorrer do processo de 
moagem, as partículas dos pós são repetidamente levadas à solda, quebra e ressolda num 
intenso processo cíclico de transferência de energia que possibilita a nanoestruturação dos 
materiais moídos [39]. 

No processo de quebra e solda das partículas certa quantidade de pó é pressionada 
pelo choque, Figura 3-1. 


esfera 


esfera 


Figura 3-1: Ilustração de uma colisão esfera-pó-esfera durante um processo de moagem de alta 
energia,adaptado da referência [39]. 


Como pode se observar, os impactos geram deformações plásticas e trituração. 
Desta forma, se duas ou mais partículas são deformadas e sobrepostas, podem se agregar 
por um mecanismo de solda a frio e assim aumentar o seu tamanho até haver uma nova fase 
de quebra. Nos estágios iniciais da moagem, as partículas possuem dureza relativamente 
baixa, existindo uma grande tendência de se reunirem e formarem aglomerados. Neste 
estágio as partículas possuem uma estrutura morfológica com várias camadas constituídas 
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de várias combinações dos elementos constituintes. 

Com a contínua deformação, as partículas ficam mais rígidas e se fraturam devido 
ao mecanismo de fadiga. Fragmentos gerados por esse processo tendem a diminuir de 
tamanho na ausência de forças de aglomeração. Até este momento, a tendência à fratura 
predomina sobre a soldagem a frio. Devido aos contínuos impactos das esferas, a estrutura 
das partículas é rapidamente refinada, porém, o tamanho das partículas continua o mesmo. 

Depois de moído por certo período de tempo, um estágio estacionário é atingido 
após alcançar um equilíbrio entre a taxa de soldagem, a qual tende a aumentar o tamanho 
das partículas, e a taxa de fratura, a qual tende a diminuir o tamanho das partículas. 
Pequenas partículas são capazes de suportar deformações sem sofrer fratura, tendem a se 
unir em partículas maiores, levando a formar partículas de tamanho intermediário depois de 
certo período de tempo. Neste estágio, cada partícula possui todos os componentes iniciais 
na proporção em que foram misturados, tonando-se homogêneas, e atingem seu limite de 
dureza devido ao acúmulo de deformações [39]. 

O processo de moagem possui vários parâmetros que influenciam nas propriedades 
do produto final. Tais como: razão massa das esferas/massa dos pós, carga do vaso de 
moagem, velocidade de rotação ou frequência de vibração, tempo de moagem, atmosfera 
de moagem e os tipos de aparelhos utilizados na moagem. 


21 PROCESSO DE SINTERIZAÇÃO 


O processo de sinterização pode ser definido como a união físico-química de um 
sistema de partículas empacotadas pela formação de ligações entre estas partículas, 
acompanhada por um fluxo de massa (Figura 3-2). Asinterização é geralmente acompanhada 
por um aumento na resistência mecânica, na ductilidade e pelo aumento da densidade [40], 
[41]. 


Fluxo de massa 


e 


p 


EA 


Figura 3-2: Modelo de esfera para duas partículas de raio r, sendo x o raio do pescoço e p o raio de 
curvatura da superfície na região do pescoço [42]. 
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A sinterização ocorre quando as partículas ficam em estreito contato e a temperatura 
é suficientemente alta para produzir a união por coalescência. Em alguns casos, quando se 
trabalha com um sistema multicomponentes e a temperatura supera a temperatura de fusão 
de algum dos componentes, tem-se a formação de uma fase líquida. A sinterização de um 
material provoca usualmente muitas mudanças nas suas propriedades. Nas cerâmicas, o 
processo de sinterização pode aumentar a resistência mecânica, a condutividade térmica, 
a densidade e possibilita a produção de peças transparentes ou translúcidas [41]. 

As mudanças nas propriedades microestruturais de um material com a sinterização 
podem ser visualizadas na Figura 3-3. 


Figura 3-3: Representação esquemática das etapas de sinterização de pós-cerâmicos, sendo (a) início 
da formação de “pescoços”, (b) fase inicial de sinterização, (c) fase intermediária de sinterização, (d) 
fase final de sinterização [43]. 


Como podemos observar na Figura 3-3(a), em temperatura adequada para que 
o fenômeno ocorra, um sistema de partículas em contato começa a formar “pescoços” 
entre si. A Figura 3-3(b) representa a fase inicial da sinterização. Na fase intermediária, 
Figura 3-3(c), a área de contato entre as partículas aumenta e os poros começam a ser 
suavizados. Não existem poros fechados nesta etapa. Na fase final da sinterização Figura 
3-3(d) os poros são fechados (não existindo comunicação entre eles) e estão localizados 
nos contornos de grão. A densificação completa ocorre quando todos os poros são fechados. 

O processo de sinterização convencional trata-se de um processo não- linear e tem 
sua fase final caracterizada pelo crescimento de grãos e surgimento de poros fechados e 
isolados nos contornos de grão. Assim, a estrutura do corpo cerâmico é vista como grãos 
contendo poros isolados nos seus contornos e não mais como partículas isoladas. Os 
poros isolados no corpo cerâmico agem como uma segunda fase impedindo o crescimento 
muito rápido dos grãos. Com o crescimento dos grãos os poros fechados são pressionados 
e há uma tendência do gás aprisionado nos poros se difundir pelo contorno dos grãos. Com 
essa difusão os poros diminuem de tamanho e a densidade real do material se aproxima 
da densidade teórica do mesmo. Em alguns materiais cerâmicos ocorre um crescimento 
exagerado de grão, sendo que, muitas vezes, o grão cresce extremamente rápido, 


englobando os poros e assim aprisionando-os dentro do grão. Poros deste tipo, no interior 


1 Região de ligação entre as partículas. 
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de grãos, não serão fechados, pois os processos difusionais são extremamente lentos e, 
consequentemente, a peça não atingirá um valor próximo à teórica [1]. 

Ao se trabalhar com sistemas multicomponentes, pode acontecer o aparecimento de 
fases líquidas durante a sinterização. Estas fases são capazes de acelerar e promover a 
densificação dos corpos cerâmicos, podendo transportar grande quantidade de material em 
menor tempo, se comparado a processos de difusão em escala atômica [1]. 


2.1 Método de Sinterização 
Para este trabalho foram utilizados dois métodos de sinterização: 
* | Sinterização convencional em atmosfera livre 
* - Sinterização rápida seguida de choque térmico 


A seguir é apresentada uma breve descrição deste método. 


2.2 Sinterização Convencional em atmosfera livre 


No procedimento de sinterização convencional a amostra é submetida a uma 
variação controlada de temperatura com taxas de alguns graus célsius por minuto. Tal 
sinterização consiste basicamente de uma rampa de aquecimento a taxa fixa, por exemplo 
5ºC/min, até atingir um patamar na temperatura de sinterização desejada e, finalmente, 
uma rampa de resfriamento, essa última com taxa de controle ou obedecendo a inércia 
térmica do forno utilizado. Neste trabalho utilizamos o resfriamento obedecendo a inercia 
térmica do forno, geralmente utiliza-se fornos resistivos para este tipo de sinterização, com 
controle ou não de atmosfera. 


2.3 Sinterização Rápida Seguida de Choque Térmico 


O processo de sinterização rápida seguida de choque térmico consiste em submeter 
a amostra a uma sinterização com altas taxas de aquecimento, por exemplo, neste trabalho 
utilizamos uma taxa de aproximadamente 20º C/s, e em seguida submetê-la a um rápido 
resfriamento, passando da temperatura de sinterização a temperatura ambiente em poucos 
segundos. 

Como exemplo de arranjo experimental empregado em processos de sinterização 
rápida seguida de choque térmico, podemos citar o empregado neste trabalho. Este 
se utilizou de um forno resistivo usualmente utilizado em processos de sinterização 
convencional, Figura 3-4 (1), porém adaptado para a realização da sinterização rápida. Na 
entrada do forno observa-se a adaptação de uma segunda porta, Figura 3-4 (3), com uma 
entrada central para a inserção da amostra para realização do processo de sinterização, 
sendo disposta sobre uma gaveta de concreto refratário, Figura 3-4 (4). A temperatura e 
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taxa de aquecimento da amostra pode ser acompanhada pela leitura do termopar. 


Termopar 
N (3) 


cad 


Figura 3-4: Arranjo experimental utilizado para sinterização rápida seguida de choque térmico. Onde 
(1) trata-se de um forno resistivo convencional controlado eletronicamente, (2) multímetro para 
acompanhamento da temperatura da amostra, (3) porta adaptada na entrada do forno especialmente 
para realização da sinterização rápida e (4) a gaveta sobre a qual é disposta a amostra durante a 
sinterização, adaptado de [16]. 


Neste processo a amostra é inserida rapidamente no forno, já aquecido, e atinge 
a temperatura de sinterização em poucos segundos. Depois de transcorrido o tempo de 
sinterização da amostra, essa é retirada rapidamente do interior do forno e transferida da 
gaveta para uma superfície a temperatura ambiente, sofrendo assim um choque térmico 
[16]. 


31 DIFRAÇÃO DE RAIOS X 


A difração de raios X constitui-se como uma importante técnica de análise estrutural 
de materiais cristalinos e tem como base o espalhamento de um feixe raios X pelos átomos 
que constituem a rede cristalina do material analisado, permitindo sua identificação e 
caracterização [6]. A difração de raios X em cristais tem como princípio que estes são 
formados por átomos ordenados de forma periódica em três dimensões, formando assim 
uma rede cristalina. Podemos descrever a rede cristalina em termos de parâmetros de 


rede, que são os comprimentos e ângulos que definem uma cela unitária, que por vez 
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gera a rede cristalina através de operações de simetria [6]. Considerando-se dois ou mais 
planos de uma estrutura cristalina, as condições para que ocorra a difração de uma onda 
(interferência construtiva ou destrutiva) dependem da diferença de caminho percorrido e 
do comprimento de onda da radiação incidente. Na Figura 3-5 são mostrados os feixes 
incidentes e espalhados. 


Figura 3-5: Esquema ilustrativo do fenômeno de difração de raios X (Lei de Bragg), adaptado da 
referência [6]. 


Quando o comprimento de onda de uma onda eletromagnética é comensurável com 
as dimensões atômicas da rede, esta onda interage com o cristal e é difratada de acordo 
com lei de Bragg [45], dada por: 


2dsinô0 =nA N= 1,23, co (2) 


Onde, À é o comprimento de onda da radiação incidente, 9 corresponde ao ângulo 
medido entre o feixe incidente e determinados planos do cristal, d é a distância entre 
os planos de átomos e n é um número inteiro (ordem de difração). Debye e Scherrer 
desenvolveram uma técnica para obter os padrões de difração de raios X para amostras 
em pó, esta técnica consiste em fazer incidir um feixe monocromático de raios X sobre a 
amostra, que possui cristais orientados de forma aleatória [6]. As direções para as quais 
ocorre difração são determinadas pela simetria da rede cristalina, que ao formar planos de 
átomos em distâncias características irão difratar essa radiação incidente em determinados 
ângulos com diferentes intensidades, gerando um difratograma padrão para cada simetria, 
como o exemplo da Figura 3-6. 
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Intensidade (u. a.) 


0 20 40 eo so 100 120 140 160 
26 (graus) 


Figura 3-6: Resultado de uma medida de difração. 


41 REFINAMENTO ESTRUTURAL: MÉTODO DE RIETVELD 


Desenvolvido por Rietveld no fim da década de 60, o refinamento estrutural para um 
difratograma de raios X começou a ser utilizado quase uma década depois [45]. O método 
para o refinamento tem como base a construção de um padrão de difração calculado a 
partir de informações fornecidas pelo usuário. O ajuste desse padrão de difração calculado 
é feito refinando simultaneamente os parâmetros instrumentais e estruturais, relacionando- 
os com as características físicas da amostra analisada, visando obter uma mínima diferença 
em comparação com o padrão de difração experimental. 

Por meio das informações cristalográficas do grupo espacial, parâmetros de rede 
e posições atômicas com valores próximos aos valores reais do material em estudo, um 
padrão de difração pode ser simulado com o uso de uma equação ou modelo, que fornece 
a intensidade de cada reflexão sugerida por Rietveld [46] na forma: 


Ya =5 D LelFel? O (28: — 200) PRA + Yi (3) 
k 


Sendo: 

s é o fator de escala 

k representa os índices de Miller, h k I, para a reflexão de Bragg. 

L, este fator contém os fatores de Lorentz, polarização e multiplicidade. 
F, é o fator de estrutura da k-ésima reflexão de Bragg. 

O é a função perfil de reflexão. 

20, é o ângulo do i-ésimo ponto no padrão. 
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20, é o ângulo de Bragg calculado. 

P, é a função orientação preferencial. 

Aéo fator de absorção. 

Y,: é a intensidade da linha de base no i-ésimo passo 


A equação 3 dá a intensidade simulada para definir um padrão de difração simulado 
que possa se ajustar da melhor forma a um padrão de difração experimental. A posição de 
cada reflexão é obtida pelos parâmetros de rede, grupo espacial e por meio da lei de Bragg. 

A aproximação entre os padrões calculado e experimental é feita pelo método dos 
mínimos quadrados. A quantidade minimizada [46] é 


S= DO — Ya)? (4) 


Ee 
Onde: w, = Yi 
y, é a intensidade experimental da i-ésima reflexão, 
7, é a i-ésima intensidade calculada. 


Deste modo, quando este resíduo for mínimo encontra-se o padrão simulado que 
melhor se ajusta aos pontos do padrão observado experimentalmente [46]. 

Para se obter um bom refinamento é necessário ter medidas precisas do padrão de 
difração além da experiência e capacidade para selecionar corretamente uma sequência de 
parâmetros a serem refinados, de maneira a se obter um modelo estrutural inicial próximo 


da estrutura real do material em estudos. 


4.4 Fator de estrutura 


O fator de estrutura F, depende da posição do átomo na célula unitária e de seu 
poder de espalhamento, ele fornece a intensidade dos picos de difração e onde os picos 
devem aparecer, este fator é dado na forma: 


Fe= : N;fiexp|2m i(ha; + dy + Iz;)| exp[-M;] (5) 
j 


Onde: 

f, é o fator de espalhamento para o átomo 

ii h k | são indices de Miller. 

X, Y; e z; são coordenadas do j-ésimo átomo da celular unitária. 
M, é o parâmetro de vibração térmico. 

Nj é a ocupação de sítio dividido pela simetria máxima do sítio. 
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4.5 Fator de escala 


O fator de escala s é a constante que ajusta a intensidade em relação à altura dos 
picos, do padrão de difração teórico e do padrão de difração experimental. A determinação 
deste fator é importante para a análise quantitativadas fases. 


4.6 Fator Lorentz, polarização e multiplicidade - L, 


O fator de Lorentz, dado pela equação 6, 
ne | 
a / 4sen20cos6 (6) 


surge da divergência dos raios X, este fator pode ser agrupado ao fator de polarização 
dado pela equação 7, 


Fp = 1+ cos*26/, (7) 


que representa a intensidade da onda espalhada por um único elétron, o efeito 
global destes fatores é descrever a intensidade das reflexões que ocorrem em ângulos 
intermediários. O fator de multiplicidade p contribui para o cálculo da intensidade como 
sendo a proporção relativa de planos contribuintes para a reflexão. O fator L, é dado pela 
equação 8: 
1+ cos226 


=p-——— 8 
Le P'senZ8coso (8) 


4.7 Função de perfil de reflexão O 


A adaptação de uma função matemática que defina adequadamente o formato de 
pico, largura e posições das reflexões de Bragg é uma etapa fundamental na obtenção do 
padrão calculado. A partir do ajuste do perfil das reflexões de Bragg são obtidas informações 
importantes contidas em um padrão de difração, por exemplo, intensidade de picos de 
reflexão, posições de picos, largura, forma do decaimento das curvas dos perfis de reflexão 
e a área integrada [46]. Entre as funções mais utilizadas nos diferentes programas para o 
refinamento Rietveld, estão as funções de perfil: pseudo-Voigt, Pearson VIl e pseudo-Voigt 
de Toraya. 

A função de perfil de reflexão leva em consideração efeitos das características 
instrumentais e possíveis características da amostra que causam aberrações como: 
absorção (transparência), deslocamento da amostra, alargamento dos perfis de reflexão 
devido ao tamanho dos cristalitos. 
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As larguras dos picos são modeladas no refinamento como uma função da tg6 e 
relaciona a largura à meia altura (H), FWHM (full width at half maximum), com a ângulo 6 
de difração, como podemos ver na equação 9. 


Hº = Utg?0+Vtgo+W (9) 
Onde U,V e W são parâmetros refináveis. 
4.8 Função Voigt 
Para se obter a forma dos picos de difração de raios X, vários tipos de funções 


foram formuladas e testadas. A função Voigt que é uma convolução das funções Gaussiana 
(equação 10) [46], 


Co Co(246;n)? 
= . pi ai 10 
G E exp ( HZ (10) 


onde: 

C, é a constante de normalização (= 4ln2); 

H, é a largura total a meia altura do pico h = (h k 1); 

240, é a distância em graus (20) da posição do pico de Bragg até o iésimo ponto 
onde a intensidade está sendo calculada. 

e Lorentziana (equação11) [46], 


EO Ci(2ABn | 
a Ss. 11 
e A, 1+ AR (11) 


onde: 

H, é a largura total a meia altura do pico h = (h k 1); 

C, é a constante de normalização (= 4); 

240,, é a distância em graus (20) da posição do pico de Bragg até o iésimo ponto 
onde a intensidade está sendo calculada. 

mostra bons resultados ao calcular tamanho de partículas e efeitos de strain. 


4.9 Função pseudo-Voigt 

A função pseudo-Voigt, PV é uma aproximação analítica da função de Voigt, esta 
função é muito recomendada devido à estabilidade de refinamento e por permitir identificar 
uma provável presença de uma distribuição bimodal de tamanho do cristalito quando n é 
maior do que um [46], [47]. Esta função é descrita por: 
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PV=nL+(1-n)G (12) 


onde mn, Le G são, respectivamente, a fração lorentziana, a função de Lorentz ea 
função de Gauss. 
A fração lorentziana é definida por, 


n=NA+NB+20 (13) 


onde NA e NB são parâmetros refináveis e devem ser positivos. 


4.10 Split pseudo-Voigt - Função de Toraya. 


A função pseudo-Voigt mostra ótimos resultados para a difração de raios X 
convencional, radiação de síncrotron e nêutrons, porém variações sutis como as que 
ocorrem na base dos picos não são bem descritas, prejudicando o ajuste global. A função 
split pseudo-Voigt modificada por Toraya possui algumas vantagens sobre a função pseudo- 
Voigt e as demais funções de perfil, ela consegue modelar algumas formas mais complexas 
de picos, dando um melhor ajuste e uma análise mais precisa, ainda pode se ajustar aos 
mais variados tipos de radiação sem alterar o modelo matemático [48], [49]. 


4.11 Função orientação preferencial P, 


A orientação preferencial surge quando os cristais da amostra tendem a estarem 
mais orientados em uma direção do que em outras [46]. A orientação preferencial altera 
as intensidades relativas das reflexões de certos planos h k | prejudicando o refinamento. 
As distorções causadas podem ser modeladas matematicamente como uma “função 
orientação preferencial”, P, que é dada por: 


Py = (Gicos?a + (1/c,) senta) fa (14) 


Onde G, é um parâmetro refinável e a é o ângulo entre dk* e o eixo de difração. 


4.12 Intensidade da radiação de fundo (Background) (y,,) 


A radiação de fundo é consequência de vários fatores, como fluorescência da 
amostra, ruído no detector, fases amorfas na amostra, espalhamento de raios X no porta 
amostra [46]. A intensidade do background pode ser obtida de uma interpolação linear de 
pontos selecionados fora dos picos no padrão. Este método é adequado quando o padrão 
é relativamente simples. Para padrões mais complexos deve ser utilizada uma função de 
background especifica. Uma função simples é um polinômio de quinta ordem provido com 
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um operador específico para permitir uma flexibilidade no ajuste de picos largos e nas 
curvas de background [46]. O modelamento é descrito pela equação 15: 


5 
20; do 
Ybi = >. Bm I( /BKPOS) = 1] (15) 
m=0 
Onde B, representa os coeficientes do polinômio, que são parâmetros ajustáveis; 
BKPOS é o operador específico usado no arquivo de controle. 


4.13 Indicadores de qualidade do refinamento 


O refinamento chega ao fim quando os parâmetros não variam mais e a função de 
minimização alcança o valor mínimo. No entanto a convergência deve ser acompanhada 
através de alguns fatores que são calculados ao final de cada ciclo, e fornece ao usuário 
subsídio para tomar decisões sobre dar continuidade, parar ou finalizar o refinamento 
[50]. Podemos considerar que o melhor resultado é aquele que fornece o difratograma de 
raios X calculado mais próximo possível do observado. Através de critérios numéricos e 
gráficos podemos confirmar quantitativamente esta aproximação, os critérios numéricos 
são conhecidos como resíduos - (fatores de qualidade de ajuste). 

Os resíduos mais usuais (em %) são: 


* Resíduo Rp 


Daly; — Yeil 


R, = 100 
á Dilyil 


(16) 


Este resíduo quantifica a diferença entre as intensidades calculadas e observadas 
experimentalmente. 


* Resíduo Rwp 


Di Wi(Yi — Yci)? (17) 


Ryp = 100 
cc Diwi yf 


onde w, é o peso atribuído ao i-ésimo ponto durante o refinamento (w, = 1,) 


O fator de qualidade de ajuste Rwp é outro indicador, que considera o resíduo em seu 
numerador e da maior peso a picos de maior intensidade. Os fatores que modificam Rwp 
são as diferenças na forma dos picos, por exemplo a largura, e a estatística da radiação 
de background. Para uma análise satisfatória dos resíduos Rp e Rwp são considerados 
valores abaixo de 10%, devido a presença de ruídos que impedem um ajuste perfeito entre 
o padrão observado e o calculado [50]. 

Outro índice importante para o análise de qualidade do refinamento é o valor 
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estatisticamente esperado para o Rwp (equação 18). 


* Resíduo Re 


W=Pr+Op” 


E WiyZ o 


R. = 100 


Onde P é o número de parâmetros refinados, N é o número de pontos experimentais 
do difratograma e € é o número de vínculos aplicados. 

O resíduo R, está diretamente relacionado com a qualidade dos dados coletados 
no difratograma, pois a presença de maior intensidade de ruído, ou picos com baixa 
intensidade eleva seu valor [50]. 

Teoricamente Re corresponde ao menor valor de Ryo! e quanto mais próximo de 
1arazão R,/R, estiver, melhor será o refinamento. A razão entre R,, e R, fornece outro 
indicador, o Goodness-of-fit ou qualidade do ajuste (GOF) ou x?, que em um refinamento 
perfeito tem valor 1. 


PA 


(GOP) = 2 = (52) (19) 


Para obtermos um refinamento com qualidade considerada satisfatória é necessário 
que x? esteja abaixo de 4, entretanto este não deve ser inferior a 1, pois valores de 7? 
menores que 1 podem indicar que R teve seu valor subestimado, ou que o refinamento 


está ajustando o ruído, o que não é desejável [50]. 


51 FULLPROF 


O programa FullProf suite (para Windows e Linux) é utilizado para fazer o refinamento 
Rietveld, através de padrões gerados por difração de nêutrons ou raios X, coletados com 
o passo de contagem constante em 20. Escrito na linguagem Fortran, ele é organizado de 
modo a se adequar em diferentes computadores. 

Entre as vantagens do programa, podemos destacar a possibilidade de escolha de 
uma função de perfil (Gaussiana, Lorentziana, Lorentzianas Modificadas 1 e 2, Pseudo- 
Voigt, Pearson-VIl ou Thompson-Cox-Hastings) diferente para cada fase e a facilidade 
quanto a visualização de gráficos (através do programa WinplotR). Outra vantagem que 
podemos destacar é a gratuidade do programa. 


61 RIETAN-FP. 


Desenvolvido por Fujio Izumi [51], o RIETAN-FP é um programa para a realização 
de refinamentos estruturais pelo método de Rietveld. Assim como o FullProf o programa 
é codificado na linguagem Fortran. Pode ser utilizado para difração de raios X, nêutrons e 
simulações. Os três sistemas principais do Rietan-FP são: 
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* Padrão de composição pelo método de Le Bail. 
* | Refinamento dos parâmetros de rede e estruturais pelo método deRietveld. 
* | Cálculo da densidade eletrônica pelo método da máxima da entropia (MEM). 


Avantagem do RIETAN-FP sobre outros programas de análise de Rietveld é a opção 
de uma nova função de perfil pseudo-Voigt modificada por Toraya e sua proximidade com 
o programa VESTA, um poderoso sistema de visualização 3D. Ainda podemos citar a sua 
estabilidade, a convergência automática para um mínimo, a velocidade de refinamento e a 
simplicidade para o manuseio do programa [51]. 


71 MÉTODO DA MÁXIMA ENTROPIA (MEM) NA CRISTALOGRAFIA 


O método da máxima entropia (MEM) é utilizado para mapear em alta resolução as 
densidades eletrônicas no interior da célula unitária. A densidade eletrônica é descrita pelo 
fator de estrutura obtido com o cálculo de difração de raios X e pelo refinamento estrutural. 

O fator de estrutura é uma função de distribuição de probabilidade eletrônica que 
descreve as posições atômicas, como podemos ver na equação 20. 


1 ! 1 
F(h,k,D = [ [ [ V p(x,y,z)exp|-2mi(hx + ky + Iz)]dxdydz (20) 
00-20 


Onde 

V é o volume da célula unitária 

p(x, y, z) é a densidade eletrônica no interior da célula unitária 

(Ah, k, 1) são os índices de Miller 

(x, y, Z) são as coordenadas atômicas dos átomos pertencentes a célula unitária 


Com a transformada inversa de Fourier obtemos a densidade eletrônica p(x, y, z), 
descrita na equação 21. 


1 
p(x,y,7) = 7 ) 22 F(hkDexp[-2mi(hx + ky + I2)] (21) 


Porém no cálculo da densidade eletrônica (equação 21), informações são perdidas, 
devido ao problema de fases dos raios X e do cálculo da transformada de Fourier, entre 
outros [52]. A partir da técnica do MEM esta informação faltante é corrigida, a técnica 
age como uma lente virtual que ajusta a imagem com defeitos por meio de métodos 
computacionais tornando-a mais limpa. Portanto o método de máxima entropia recupera 
a informação que está faltando permitindo a reconstrução da densidade eletrônica de alta 
resolução em imagem tridimensional. 

A aplicação do MEM associado com o método de Rietveld em dados cristalográficos, 
permite maior precisão em relação a posição atômica e melhor visualização de possíveis 
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deformações causadas por dopagem de amostra. 

O princípio geral da análise do MEM é encontrar o máximo de informação relativa 
à entropia, S, sob algumas restrições de procedimentos iterativos [53]. A densidade de 
alta resolução é calculada pela otimização da função de entropia de informação dada pela 
equação 22. 


S = —D pulnde (22) 


Onde p, é densidade normalizada na posição r,, e r, é a densidade normalizada 
derivada de informações prévias para a posição r,. 


81 DYSNOMIA 


O programa Dysnomia é utilizado para fazer análise do MEM, este programa é 
sucessor do PRIMA e faz parte do sistema de visualização 3D VENUS. Ele é escrito na 
linguagem de programação C++ por meio de avanços na análise do MEM para monocristais 
e dados de difração de pó [53]. 

O Dysnomia oferece maior velocidade requerendo menos memória em cálculos, 
tornando-o adequado para a execução em computadores pessoais. 


91 LIGAÇÕES QUÍMICAS 


Geralmente os sistemas químicos tendem a um estado de menor energia, de forma 
natural, assim podemos dizer que os átomos se unem para ocupar um estado de menor 
energia, pois quando ligados, estes possuem menor energia do que quando separados. Os 
mapas de densidade eletrônica juntamente com o cálculo da mínima densidade eletrônica 
(MED) permitem avaliarmos as ligações químicas entre os átomos. Entre as ligações 
químicas podemos destacar a iônica e a covalente, alvo de nossos estudos. 

Átomos e íons possuem uma densa nuvem eletrônica, onde ao analisarmos 
as ligações podemos perceber se há ou não sobreposição das nuvens eletrônicas. 
A sobreposição caracteriza uma ligação covalente enquanto que a não sobreposição 
evidencia uma ligação iônica. 


10 | LIGAÇÃO IÔNICA 

Os átomos ao perder (ou ganhar) elétrons se tornam íons com carga positiva (ou 
negativa), chamado de cátion (ou ânion) respectivamente. 

Uma vez que os íons estão formados, a força eletrostática entra em ação fazendo 
os cátions atraírem os ânions, e vice-versa. Para aumentar ao máximo as atrações cátion- 
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ânion, e diminuir tanto quanto possível as repulsões cátion- cátion e ânion-ânion, cada íon 
positivo fica rodeado pelo maior número possível de íons negativos, por sua vez os íons 
negativos fazem o inverso, sendo envolvidos por íons positivos, como podemos ver na 
Figura 3-7. 


c 


Figura 3-7: Disposição dos íons em uma célula unitária tipo perovskita. 


Em nosso trabalho ao analisarmos as distâncias entre os íons consideramos a 
tabela? de Shannon [54] sobre os raios iônicos. Para reproduzir distâncias interatômicas 
Shannon levou em consideração o número de coordenação, a rotação eletrônica, a 
covalência, as forças de repulsão e a distorção poliédrica. Com estas limitações, os raios 
são independentes do tipo de estrutura. Tanto para o cátion como para o ânion os raios 
variam de acordo com número de coordenação [54]. 


10.14 Ligação Covalente 


A ligação covalente é caracterizada pela partilha de elétron entre os átomos, ou 
seja, a sobreposição de nuvens eletrônicas, indicando uma ligação de covalência. O 
composto de (Bi, Nd,)(Fe,,Co)JO,, apesar de formado por íons, apresenta ligações de 
caráter covalente, este comportamento tem sido apontado como propulsor do surgimento 
da polarização elétrica [48]. 

Os Lone Pairs (pares de elétrons não compartilhados), presentes no bismuto 
(Bi) são mediadores em ligações covalentes e possivelmente causadores de distorções 
estruturais. Na origem do ferroismo em compostos multiferróicos as ligações de caráter 
covalente e os Lone Pairs existentes desempenham papel fundamental na ordenação 


2 Atabela com os raios iônicos e número de coordenação encontra-se no anexo |. 
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elástica e magnética [55]. 

No presente trabalho analisamos as densidades eletrônicas, obtidas pelo método 
de máxima entropia, entre os íons do composto (Bi, Nd,)(Fe, Co JO, visando entender o 
comportamento das ligações químicas entre os íons de (Bife O )e (Fe te 0). 
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EXPERIMENTAL 


Neste capítulo serão apresentadas as técnicas para o processamento do composto 
(Bi, Nd (Fe, ,Co JO, para os valores de (x = 0,0;0,1 e 0,2) e (y = 0,0; 0,01 e 0,02), tais 
como a preparação das amostras, difração de raios X, refinamento estrutural, cálculos de 
densidade eletrônica e de mínima densidade eletrônica. 


11 PREPARAÇÃO DAS AMOSTRAS 


No processo de preparação das amostras do composto (Bi, Nd)(Fe,,CojO, 
foram utilizados os precursores Fe,O, (Aldrich, pureza analítica de 99, 98%), (Bi,O, Vetec, 
pureza analítica de 99, 8%), (Nd,O, Aldrich, pureza analítica > 99,9%) e (Co,O, Aldrich, 
pureza analítica > 70%). Todos os pós precursores foram pesados em balança analítica e 
misturados em proporções de acordo com a estequiometria desejada, as composições do 
composto (Bi, Nd)(Fe,,Co)O, escolhidas foram (x=0,0;0,1e0,2)e (y=0,0;0,01 e 0,02). 
O vaso de moagem utilizado, assim como as esferas, foram de aço/carbono endurecido 
V C-131.A moagem em altas energias foi realizada num moinho planetário Retsch PM 
100, em atmosfera ambiente seca. A razão massa das bolas/massa dos óxidos foi de 20:1, 
sendo 2 g a massa dos óxidos, a velocidade de moagem foi de 250 rpm por uma hora. Para 


todas as amostras foram utilizadas as mesmas condições. 


21 TRATAMENTOS TÉRMICOS 


Após a moagem de alta energia do composto de (Bi, Ndj(Fe,,CoJO, os pós 
obtidos receberam o tratamento térmico adequado para eliminar fases espúrias, para isso 
utilizamos dois tipos de sinterização: a convencional em atmosfera livre (C) e a rápida 
seguida de choque térmico (RSCT). 

Na Tabela 4-1 vemos os tipos de tratamento térmico que cada amostra recebeu para 
se obter a formação de fases. 


Dopagem Tipo de Sinterização Temperatura (ºC) Tempo (Min) 
x=0,0ey=0,01 RSCT 870 E] 
x=0,0ey=0,02 RSCT 880 3 
x=0,1ey=0,0 (0) 830 60 
x=0,22ey=0,0 (6, 850 60 
x=0,1ey=0,01 RSCT 900 4 
x=0,2ey=0,02 RSCT 900 6 


Tabela 4-1: Tratamentos térmicos para as dopagens do composto (Bi, Nd (Fe, ,Co)JO.. 
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Os passos seguidos para as sinterizações foram previamente reportados na 
referência [1]. 


31 DIFRATOMETRIA DE RAIOS X 


A difratometria de raios X do composto (Bi, Nd (Fe, Co, JO, foi realizada utilizando 
o difratômetro Shimadzu modelo XRD-7000, com radiação de Cu Ka do laboratório do 
grupo de desenvolvimento de dispositivos multifuncionais GDDM — UEM. As medidas 
foram realizadas o passo de 2º/min, de 10º até 130º, em modo contínuo, todas realizadas à 
temperatura ambiente. Os difratogramas obtidos por difração de raios X foram utilizados na 
identificação de fases formadas durante o processo de sinterização. Desta forma tivemos 
que repetir o processo de sinterização diversas vezes até conseguir as fases desejadas. 


41 REFINAMENTO ESTRUTURAL E CÁLCULO DE DENSIDADE ELETRÔNICA 

Os padrões de difração foram analisados por refinamento estrutural, método de 
Rietveld, utilizando os programas FullProf e Rietan-FP, o motivo do uso dos dois programas 
está descrito no capítulo 5. As densidades eletrônicas foram calculadas pelo programa 
Dysnomia e visualizadas no programa VESTA. 
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RESULTADOS 


Ao longo deste capítulo, apresentaremos os resultados obtidos sobre as estruturas 
cristalinas de composições do sistema multiferróico (Bi, Nd)(Fe,,Co)O, através do 
refinamento estrutural (método de Rietveld) e por cálculos de densidades eletrônicas 
realizados pelo Método da Máxima Entropia (MEM). 


11 REFINAMENTO ESTRUTURAL DE RIETVELD 


Nesta sessão são apresentadas as análises dos padrões de difração de raios X 
do sistema multiferróico (Bi, Nd)(Fe,,Co)O, (x=0,0;0,1€e0,2)e(y=0,0; 0,01 e 0,02). 
Inicialmente utilizamos o programa FullProf para o refinamento de Rietveld, fazendo uso 
da função Pseudo-Voigt e Background com interpolação de pontos. Como sugerido na 
literatura, foram utilizados os grupos espaciais Romboédrico (R3c) e Ortorrômbico (Pbnm) 
(fases individuais) e o modelo de coexistência das duas fases para o refinamento. O fator 
de qualidade X obtidos para os refinamentos estão listados na Tabela 5-1. 


Dopagem Grupo Espacial 
dE R3c Pbnm R3c&Pbnm 
x=0,0ey=0,01 7,8 47,2 9,9 
x=0,0ey=0,02 4,7 40,9 3,9 
x=0,1ey=0,0 5,1 19,5 4,9 
x=0,2ey=0,0 Sil 23,4 4,9 
x=0,1ey=0,01 5,2 14,4 6,4 
x=0,2ey=0,02 8,9 Td 10,5 


Tabela 5-1: Valores de qualidade de refinamento x?, obtidos com o programa FullProf. 


Mesmo atingindo um X <4 com a coexistência das duas fases para a dopagem com 
x=0,0ey= 0,02, 0 programa FullProf ao realizar o refinamento utilizou aproximadamente 
100% da fase R3c para o ajuste, o que não caracteriza uma transição estrutural. Como 
podemos, não foi obtida uma boa qualidade nos refinamentos, devido a esse fato optamos 
por refazer os refinamentos utilizando o programa Rietan-FP. 

Com o programa Rietan-FP e a função Split Pseudo-Voigt modificada por Toraya, 
conseguimos uma melhora significativa no refinamento de Rietveld. Adotamos os mesmos 
modelos de grupos espaciais utilizados no programa FullProf, e ainda fizemos uma 
tentativa com a coexistência dos grupos espaciais Romboédrico (R3c) e Triclínico (P1) para 


o compósito (Bi, ;Nd, »)(Fe,,92006,02)05+ como veremos no final desta sessão. 


0,02) 
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1.1 Dopagem no Sítio B - Bi(Fe,,Co)O, 


Iniciamos com a análise do difratograma do composto Bi(Fe, 900,,01)0, (Figura 
5-1), para esta amostra o melhor ajuste obtido foi com o grupo espacial R3c, para melhor 
compreensão deste e dos demais difratogramas que vierem por seguinte, descreveremos a 
legenda, na qual a linha continua azul esboça o difratograma teórico, ou seja, o difratograma 
a ser ajustado pelo refinamento de maneira a ficar o mais próximo possível do difratograma 
experimental, os círculos vermelhos descrevem o difratograma obtido experimentalmente, 
os traços verticais apontam as posições de Bragg e pôr fim a linha continua cinza nos dá a 


diferença do difratograma obtido experimentalmente em relação ao teórico. 
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Figura 5-1: Difratograma de raio X para amostra de Bi(Fe, C0,01)0,» refinada pelo método de Rietveld 
com o grupo espacial R3c. 


Os resultados com os parâmetros estruturais da amostra Bi(Fe0,99C00,01)03 estão 
listados na Tabela 5-2. 


Grupo Espacial R, Rio Chi? 
R3c 5,59 10,64 3,70 
a=b=5,578 À c=13,867 À 
íon x y Z; 
Bi 0 0 0,0200 
Fe**/Co** 0 0 0,2200 
o” 0,4221 -0,0546 0,9778 
Tabela 5-2: Dados estruturais da fase Romboédrica obtidos pelo método de Rietveld para a amostra 
Bi(Fe0,99C00,01)08. 


A Qualidade do refinamento pode ser verificada pelos parâmetros de concordância, 
[o) Y obtido foi 3,70. 
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Na Figura 5-2 podemos observar uma boa concordância entre o difratograma 
experimental e o teórico da amostra Bi(Fe0,98C00,02)058. 
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Figura 5-2: Difratograma de raio X para amostra de Bi(Fe, ,:C0,02)0.» refinada pelo método de Rietveld 
com o grupo espacial R3c. 


A seguir podemos verificar a qualidade do refinamento e os parâmetros estruturais 
na Tabela 5-3. 


Grupo Espacial R, Rio Chi? 
R3c Snl2 8,23 2,58 
a=b=5,578Á c=13,866 À 
fon x y F4 
Bi 0 0 -0,0002 
Fe*:/Co* 0 0 0,2203 
o” 0,4374 -0,0232 0,9627 


Tabela 5-3: Dados estruturais da fase Romboédrica obtidos pelo método de Rietveld para a amostra 
Bi(Fe, 800, o2)O,: 


0,02 

Assim como na amostra anterior a estrutura que melhor representa o composto 
Bi(Fe, 9s€06,02)0, é a estrutura romboédrica com grupo espacial R3c. O resultado estrutural 
para o BiFeO, dopado com Co, obtido em nosso trabalho, segue o mesmo da referência 
[37], na qual descreve que a estrutura do composto dopado é muito semelhante à do 
BiFeO, puro, e que possui estrutura romboédrica (R3c). 
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1.2 Dopagem no Sítio A - (Bi, Nd,)FeoO, 
Segundo as referências [23], [32]-[36] listadas na Tabela 2-3 quando feito a dopagem 
com x < 0,2, a estrutura mantém-se romboédrica, o resultado encontrado neste trabalho 


concorda com o proposto pelos autores. A Figura 5-3 esboça o difratograma de raio X para 
amostra de (BiyoNd,,)FeO,, refinada com o grupo espacial R3c. 
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Figura 5-3: Difratograma de raio X para amostra de (Bi, ;Nd, )FeO,, refinada pelo método de Rietveld 
com o grupo espacial R3c. 


Os resultados com os parâmetros estruturais encontram-se na Tabela 5-4, através 
dos parâmetros de concordância podemos verificar a boa qualidade do refinamento, pois 
2 
obtivemos um x igual a 3,60. 


Grupo Espacial R, Rio Chi? 
R3c 5,26 10,00 3,60 
a=b=5,571 À c=13,801 À 
íon x y F4 
Bi/Nd*+ 0 0 -0,0075 
Fe+ 0 0 0,2215 
o? 0,4395 -0,0013 0,9594 


Tabela 5-4: Dados estruturais da fase Romboédrica obtidos pelo método de Rietveld para a amostra 
(Bi,oNd,,)Feo,. 


Segundo os autores [23], [32]-[36] listados na Tabela 2-3, para a amostra de 


(Bis sNd,,)FeO, há uma grande divergência em relação a estrutura que melhor represente-a, 
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os autores da referência [33] descreveram a coexistência de fases romboédrica e 
ortorrômbica para esta amostra, os autores das referências [23], [32], [34] corroboram que 
a estrutura que melhor representa a amostra de (Bi, jNd, ,)FeO, é a estrutura com simetria 
triclínica , enquanto que os autores da referência [36] exibem uma estrutura com simetria 
ortorrômbica para esta amostra. 

Na Figura 5-5 podemos observar um bom ajuste do difratograma com o grupo 
espacial R3c. 
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Figura 5-4: Difratograma de raio X para amostra de (Bi, ,Nd,,)FeO,, refinada pelo método de Rietveld 
com o grupo espacial R3c. 


Na Tabela 5-5, estão listados os parâmetros estruturais e o parâmetro de qualidade 
de refinamento, no qual obtivemos uma boa concordância, com chi? igual a 3,94. 


Grupo Espacial R, Ryo Chi? 
RSc 5,18 10,30 3,94 
a=b=5,573 À c=13,807 À 
íon x y z 
Bi**/Nd*+ 0 0 0,0029 
Fe” 0 0 0,2247 
o” 0,4405 0,0033 0,9626 


Tabela 5-5: Dados estruturais da fase Romboédrica obtidos pelo método de Rietveld para a amostra 
(Bi,aNd, )JFeo.. 


Em nosso trabalho, obtivemos um resultado análogo ao trabalho dos autores da 
referência [35], no qual verificaram que a amostra de Bi, Nd,FeO, com x < 0,2 deve ser 
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descrita por simetria romboédrica. 


1.3 Dopagem Simultânea dos Sítios A e B - (Bi, Nd,)(Fe,,Co))O, 

De acordo com os autores da referência [38], para o composto (Bi, Nd)(Fe, CojJO, 
com (x = 0,06; 0,09;0,12e 0,15) e (y = 0,03), houve uma transição de simetria romboédrica 
para triclínica com o aumento da dopagem com neodímio. Em nosso trabalho preparamos 
duas amostras com dopagem simultânea para o composto (Bi, Na (Fe, CoJO,, uma com 
x=0,10ey=0,01 e outra com x = 0,20 e y = 0,02. Para o composto (Bi, Nd, (Fes aC0,01) 
O, obtivemos uma boa concordância no refinamento do padrão de difração de raios X 
utilizando o grupo espacial R3c. Como podemos ver na Figura 5-5. 
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Figura 5-5: Difratograma de raio X para amostra de (Bi, Nd, ,)(Feç s9C0501) Os, refinada pelo método 

de Rietveld com o grupo espacial R3c. 


Os resultados com os parâmetros estruturais estão representados na Tabela 5-6. 
Podemos confirmar a qualidade do refinamento através dos parâmetros de concordância, 
para este refinamento obtivemos um x? igual 3,55. O resultado obtido para esta amostra ficou 
dentro do esperado se compararmos com o resultado obtido pelos autores da referência [38], 
que trabalharam com valores de dopagens próximos ao nosso. 
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Grupo Espacial R R Chi? 


e wp 
R3c 5,24 9,89 3,55 
a=b=5,571 À c=13,803 À 
íon x y Ez 
Bi/Nd* [o 0 0,0020 
Fe* o) 0 0,2257 
o” 0,4357 -0,0142 0,9646 


Tabela 5-6: Dados estruturais da fase Romboédrica obtidos pelo método de Rietveld para a amostra 
(Bi, Nd, (Fes gs 0001)0a- 


Para nossa amostra com x = 0,20 e y = 0,02 podemos afirmar que houve uma 
mudança na simetria do cristal, no entanto, não foi possível identificar o grupo espacial 
que melhor represente a estrutura, foram realizados refinamentos com os grupos espaciais 
romboédrico (R3c), ortorrômbico (Pbnm) e com a coexistência das duas fases da mesma 
forma que as amostras anteriores. Como não houve sucesso, ainda fizemos o refinamento 
com o grupo espacial triclínico (P1) e com a coexistência das fases R3c e P1. 

Na Figura 5-6, podemos ver o difratograma de raios X da amostra (Bi, Nd,,) 
(Fe, 950060205» refinada com o grupo espacial R3c. 
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Figura 5-6: Difratograma de raio X para amostra de (Bi, ;Nd, »)(Fe, 98 C05,02)05» refinada pelo método de 
Rietveld com o grupo espacial R3c. 


O difratograma de raio X para amostra (Bi, Nd, ,)(Fe, 9300, 02)0 


coexistência das fases R3c e Pbnm, pode ser visto na Figura 5-7. 
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Figura 5-7: Difratograma de raio X para amostra (Bi, NA, »)(Fe, 93005020,» refinada pelo método 
de Rietveld com a coexistência das fases R3c e Pbnm, sistemas Romboédrico e Ortorrômbico, 
respectivamente. 


Na Figura 5-8, temos o resultado para o difratograma de raios X refinado com a 
coexistência das fases R3ce P1. 
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Figura 5-8: Difratograma de raio X para amostra (Bi, NA, »)(Fe, 93005,02)0,» refinada pelo método de 
Rietveld com a coexistência das fases R3ce P1, sistemas Romboédrico e Triclínico, respectivamente. 


Podemos observar que os últimos três gráficos não possuem uma boa concordância 
entre o padrão de difração teórico e o padrão de difração experimental. Isso pode ser 
comprovado analisando-se a Tabela 5-7, nela notamos que os parâmetros de concordância 
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possuem valores elevados e que o x? não ficou próximo de quatro, o que indica que o 
refinamento não tem boa qualidade. 


Grupo Espacial R, Ryo Chi? 
R3c+Pbnm 5.63 15.28 7,35 


Tabela 5-7: Parâmetros de qualidade do refinamento para o composto (Bi, ;Nd, .)(Fe, 00,02) 


21 DISTÂNCIA DE LIGAÇÃO 


A Figura 5-9 mostra uma célula unitária do tipo Perovskita do composto (Bi, ,Nd,) 
(Fe, ,CojJO, construída com os dados do refinamento estrutural. 


Figura 5-9: Célula unitária tipo Perovskita do composto (Bi, Nd, (Fe, Co))O,. 
As distâncias entre os íons Bi**/Nd** — O? e Fe*/Co** — O? para o composto 


(Bi, Nd)(Fe,,CojO, estão listados nas tabelas 5-8, 5-9, 5-10, 5-11 e 5-12 de acordo com 
as dopagens. 
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Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (A) Ligação (A) 
06 — Bil dt 3,206 Reco? de 1,868 
06 — Bi2 d2 2,588 Fe —-O3 dz 1,868 
06 — Bi3 [ok 2,563 Fe — 04 dg 2,165 
06 — Bi4 dg 3,137 Fe —-05 dg 2,165 
Fe Oi ds 1,868 Resos dio 2,165 
Tabela 5-8: Distância de ligação obtida para o refinamento estrutural da amostra Bi(Fe, s600501)03: 


Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (A) Ligação (A) 
06 — Bil di 3,209 Reno 2 d6 1,851 
06 — Bi2 do 2,559 Fe —-O3 dz 1,851 
06 — Bi3 ds 2,505 Fe — 04 dg 2,153 
06 — Bi4 da 3,117 Fe — 05 dg 2,153 
Reoi ds 1,851 Reo6 dio 2,153 
Tabela 5-9: Distância de ligação obtida para o refinamento estrutural da amostra Bi(Fe, sC0, 02)03: 


Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (À) Ligação (A) 

06 — Bil [oh] 3,206 Fe-—- 02 de 1,909 

06 — Bi2 d2 2,494 Fe— 08 dz 1,909 

06 — Bi3 ds 2,473 Fe— 04 dg 2,110 

06 — Big da 3,152 Fe—O5 dg 2,110 

Fe-—O1 ds 1,909 Fe — 06 d10 2,110 


Tabela 5-10: Distância de ligação obtida para o refinamento estrutural da amostra (Bi, Nd,,)JFeo,. 


Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (A) Ligação (À) 
06 — Bil d1 3,305 Fe-— O2 d6 1,930 
06 — Biz d2 2,508 Fe-— 03 dz 1,930 
06 — Bi3 (oh) 2,387 Fe— 04 dg 2,094 
06 — Bi4 d4 3,177 Fe— 05 dg 2,094 
Fe-01 ds 1,930 Fe — 06 d1o 2,094 


Tabela 5-11: Distância de ligação obtida para o refinamento estrutural da amostra (Bi, jNd,,)FeO,. 
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Ligação Símbolo Distância de Ligação Símbolo Distância de 
Ligação (A) Ligação (A) 
06 — Bil d7 3,224 Fe — O2 de 1,862 
06 — Bi2 d2 2,521 Fe-— 03 dz 1,862 
06 — Bi3 ds 2,470 Fe — 04 dg 2,147 
06 — Bi4 d4 3,147 Fe—O5 dg 2,147 
Fe-O1 ds 1,862 Fe — 06 d1o 2,147 


Tabela 5-12: Distância de ligação obtida para o refinamento estrutural da amostra (Bi, ;Nd,,,) 
(Fe, 900, 0, )O- 


0,9900501 
O refinamento estrutural apresentou dois tipos distintos de vizinhos mais próximos na 
ligação Fe — O, três ligações longas e três ligações curtas, podemos verificar nas tabelas 
5-8, 5-9, 5-10, 5-11 e 5-12 que as distâncias d=d=d *d, =d,=d Na tabela 5-13 verificamos 
que os parâmetros de rede diminuíram com a dopagem de Na, quando comparados com os 
parâmetros de rede do BiFeoO, (tabela 2-1), indicando uma contração na estrutura. 


Composto (B Bi,,Nd)FeO, — Romboédrico (R3c) 
Dopagem Parâmetros de Rede 
a=b (À) c(Â) 
x=0,0 5,587 13,867 
x=0,1 ESTA 13,801 
x=0,2 5,578 13,807 


Tabela 5-13: Comparação dos parâmetros de rede do composto (Bi, Nd)FeoO,. 


A dopagem com neodímio causou uma diminuição no tamanho médio dos grãos de 
cerâmica. A concentração de neodímio tem um efeito sobre o tamanho de partícula, pois o 
raio iônico do neodímio é menor que o raio iônico do bismuto. 

Na tabela 5-14 podemos comparar as distâncias de ligação Fe-O e Bi-O entre 
as dopagens do composto (Bi, Nd)(Fe,,CoJO, eo BiFeo, puro. Verificamos que nas 
dopagens com Co, as três ligações longas aumentaram e as três ligações curtas diminuíram, 
comparando-se com as distâncias apresentadas para o BiFeO,, podemos dizer que houve 
uma distorção do octaedro (FeO,). 


1 Dados retirados da referência [28] 
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Distâncias de ligação entre os íons do 
Ligação Símbolo Composto (Bi, Nd,)(Fe, Co JO, 
x=0,0ey | x=0,0ey |x=0,0ey | x=01ey |x=0,2ey | x=0,1€e 
= 0,0 = 0,01 =I0102 =0,0 E (0)(0) y=0,01 
06 -— Bi1 d7 3,400 3,206 3,209 3,206 3,305 3,224 
06 — Bi2 d2 21528 2,588 2,559 2,494 2,508 Pib2i 
06 — Bi3 ds 2,310 2,563 2,505 2,473 2,387 2,470 
06 — Bi4 d4 3,210 SIS 3,117 3,152 EITA 3,147 
Fe—O1 ds 1,940 1,868 1,851 1,909 1,930 1,862 
Fe —- 02 d6 1,940 1,868 1,851 1,909 1,930 1,862 
Fe-— 083 dz 1,940 1,868 1,851 1,909 1,930 1,862 
Fe— 04 dg 2,120 2,165 2158 2,110 2,094 2,147 
Fe—O5 dg 2,120 2,165 2,153 2,110 2,094 2,147 
Fe — 06 d10 2,120 2,165 2158 2,110 2,094 2,147 
Tabela 5-14: Distâncias de ligação entre os íons do composto (Bi, Nd,J(Fe, ,Co)O, 


Em nosso trabalho consideramos os raios iônicos como definidos por Shannon [54]. 


Assim, os raios? iônicos dos íons Bit, Nd*, Feº+, Co% e O? são 1,47 À, 1,41 À, 0,78 À, 
0,75 À e 1,26 Á, respectivamente. Desta forma, se considerarmos os íons como esferas 
rígidas, as ligações entre esses íons devem ter comprimentos como os apresentados na 
tabela 5-15. 


Distâncias de ligação entre os íons do 
Ligação Composto (Bi, Nd,(Fe,,Co)O, (À) 
x=0,0ey x=0,0ey= | x=00ey=| x=0,1ey x=0,2ey x=0,1 ey 
=0,0 0,01 0,02 =0,0 = (0/0) =0,01 
Bi—O 2,173 2,13 2,13 2,72 2,11 2,12 
Fe—O 2,04 2,04 2,04 2,04 2,04 2,04 


Tabela 5-15: Distâncias de ligação calculadas com base no raios iônicos de Shannon. 


Assim, em nosso trabalho consideraremos ligações com distâncias menores que 
as apresentadas na tabela 5-15 como sendo covalentes. Por outro lado, ligações com 
distancias maiores que as da tabela 5-15 serão consideradas iônicas. Portanto, comparando 
as tabelas 5-15 e 5-14, podemos ver que algumas ligações apresentam caráter iônico e 


2 Atabela de Shannon não faz referência sobre o raio iônico do íon de Bi, 
valor utilizado fizemos um cálculo proporcional ao raio iônico do íon de No” . 


com coordenação 12. Parachegarmos ao 
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outras caráter covalente. Ainda podemos notar que as substituições com Nd e Co alteram 
o caráter de algumas ligações. 


31 DENSIDADES ELETRÔNICAS 


Nesta sessão apresentaremos os mapas das densidades eletrônicas do composto 
(B, Nd)(Fe,,Co)O,, obtidos a partir dos resultados do refinamento estrutural, utilizando 
o programa Dysnomia, o qual realiza cálculos de densidade eletrônica pelo Método de 
Máxima Entropia (MEM). Os cálculos foram realizados com resolução de 128x128x128 
pixels. Os resultados dos cálculos foram visualizados no programa VESTA. 

A análise dos mapas de densidade eletrônica nos fornecem informações valiosas 
sobre como estão distribuídos os elétrons em uma amostra cristalina e a sua relação com 
as ligações químicas. 

Nosso objetivo é analisar como as estruturas atômicas se modificam com as 
dopagens propostas neste trabalho, também usamos a técnica da Mínima Densidade de 
Eletrônica (MED) para avaliar se uma ligação possui caráter iônico ou covalente. 


3.4 Análise dos Mapas de Densidade Eletrônica 


Selecionamos os planos (001), (002), (100) e (200) da célula unitária tipo perovskita 
(Figura 5-9) para analisarmos como as interações foram se modificando a partir das 
dopagens. Nos compostos (B,,Ndj(Fe,,Co)O, as interações de ligação ocorrem entre os 
sub-níveis dos orbitais Bi” (6s)/Nd * (41) e O” (2p) e dos orbitais Fe“/Co” (3d) e O” (2p). 
A partir dos valores da mínima densidade eletrônica (MED) classificamos o caráter das 
ligações entre os íons citados. 

Em nosso trabalho ao analisarmos os mapas de densidade eletrônica e o cálculo 
de mínima densidade eletrônica, definimos que ligações com MED < 0,1 eÂ' como sendo 
iônicas, com 0,1 CÀ"< MED < 0,2 À” como iônicas com algum caráter covalente e MED 
> 0,2 €À“ como covalentes. 

Na Figura 5-10 vemos os mapas de densidade eletrônica para a amostra 
Bi(Fe, C0,01)04: Nos planos (100) e (001) podemos observar as interações dos íons Be 


0,99 0,01 
o”, e nos planos (200) e (002) as interações dos íons Fel eo”. 
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(a) (100) (b) (200) 


Figura 5-10: Mapas de contorno bidimensionais de densidade eletrônica da amostra Bi(Fe 


CO, o1)Os: 


0,99 


Os planos que possuem as interações entre os íons de Biº* e OZ, como exemplo 
os planos (100) e (001), possuem as mesmas características de ligação. Para não sermos 
repetitivos faremos nossa discussão quanto às ligações dos íons de Biº* e O? utilizando o 
plano (001). 

Na figura 5-11 vemos variação da densidade eletrônica na linha de conexão dos íons 


adjacentes da amostra Bi(Fe,. Co 


99900, 01/04: Os valores da mínima densidade eletrônica para 


as ligações desta amostra são apresentados na Tabela 5-16, que nos permite uma análise 
mais precisa quanto às ligações visualizadas na Figura 5-10. De acordo com os valores de 
MED (tabela 5-16) obtida para as ligações do plano (001) (figura 5-10(c)), podemos dizer 
que as ligações d' e d? são iônicas, que a ligação dº é iônica com caráter covalente e para 
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a ligação d* temos os Lone Pairs (covalente). Nos planos (200) e (002) (figura 5-10(b) e (d)) 
vemos as interações entre os íons de Fe*+ e OZ, as ligações dº, de d” são covalentes, já 


as ligações dº, dº e d'º são iônicas com caráter covalente. 


“85 1.0 1.5 2.0 2.5 3.0 
d(Ã) 
Figura 5-11: Variação da densidade eletrônica na linha de conexão dos íons adjacentes na amostra 
Bi(Fe, 900, o1)O5: 
Ligação Símbolo MED (CÂ?) Ligação Símbolo MED (À?) 
06 — Bil di 0,087 Fe — O2 de 0,475 
06 — Biz d2 0,080 Fe — 03 d7 0,470 
06 — Bi3 as 0,160 Fe — 04 dg 0,154 
06 — Big dá 0,243 Fe — 05 dg 0,154 
Remo ds 0,493 Fe — 06 d1o 0,116 


Tabela 5-16: Mínima Densidade Eletrônica (MED) calculada para a amostra Bi(Fe, 500, 1)Os: 


Na Figura 5-12vemosos mapas de densidade eletrônica para a amostra Bi(Fe, 5500602) 
O,, é perceptível que o aumento na porcentagem de cobalto causou mudanças de ligação 


entre os íons. 
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(a) (100) 


(a! 


O,. 


Figura 5-12: Mapas de contorno bidimensionais de densidade eletrônica da amostra Bi(Fe, ,sC0, 02) 

Analisando as interações entre os íons de Bite o” (figura 5-12(c)) e com a ajuda 
dos valores de MED (tabela 5-17) é possível classificar a ligação d! como iônica com 
caráter covalente ou como Lone Pairs e as ligações d, d, e d, como iônicas. Nos planos 
(200) e (002) (figura 5-12(b) e (d)) é possível visualizar uma maior interação do íon de ferro 
(Fe) com os íons de oxigênio (01, 02 e 03). Esta análise pode ser confirmada verificando 
os valores de MED (tabela 5-17), na qual temos a confirmação de que todas as interações 
entre os íons de Fo e O” são covalentes, com d, d,e d mais acentuados. Na figura 
5-13 vemos variação da densidade eletrônica na linha de conexão dos íons adjacentes da 
amostra Bi(Feº*Coºº2)0º. 
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d(Á) 
Figura 5-13: Variação da densidade eletrônica na linha de conexão dos íons adjacentes na amostra 
Bi(Fe, 500, o2)0.: 

Ligação Símbolo MED(€Á') | Ligação Símbolo  MED(€Ã) 
06 — Bil di 0,155 Fe — O2 d6 0,914 
06 — Bi2 d2 0,070 Fe —-0O3 d7 0,890 
OG BIS ds 0,051 Fe — 04 dg 0,300 
06 — Bi4 dq 0,047 Fe —O5 dg 0,310 
Fe -—OÍ ds 0,914 Fe —06 d10 0,316 


Tabela 5-17: Mínima Densidade Eletrônica (MED) calculada para a amostra Bi(Fe, ,,C0, q2)O..: 
Para a amostra (BiçgNd,,)Feo, exibimos os mapas de densidade eletrônica na 
Figura 5-14, para o plano (001) (figura 5-14(c)) vemos que apesar do íon de oxigênio estar 
deslocado à esquerda, a nuvem eletrônica se mostra simétrica em relação aos íons de 
bismuto, em primeira análise podemos dizer que as ligações entre estes é iônica. Nos 
planos (200) e (002) (figura 5-14(b) e (d)) a ligação covalente entre os íons de Fete o” 


fica evidente. 
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(a) (100) (b) (200) 


(d) (002) 


e(6) 


Figura 5-14: Mapas de contorno bidimensionais de densidade eletrônica da amostra (Bi, Nd, ,)Feo,. 


Na figura 5-15 vemos variação da densidade eletrônica na linha de conexão dos 
íons adjacentes da amostra (Bi, Nd, )JFeo.. A análise da mínima densidade eletrônica para 
a amostra (Bi, Nd,,)JFeo, (Tabela 5-188) confirma queas ligações entre os íons de bismuto 
e oxigênio são iônicas extinguindo os Lone Pairs. Enquanto que as ligações entre os íons 


de ferro e oxigênio se tornaram todas covalentes com valores muito próximos entre si. 
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3.0 


Figura 5-15: Variação da densidade eletrônica na linha de conexão dos íons adjacentes na amostra 
(Bi, oNd, )Feo.. 


Ligação Símbolo MED (CÃ) Ligação Símbolo MED (CÂ) 
06 — Bi [oh] 0,005 Fe-— O2 d6 0,581 
06 -— Bi2 d2 0,003 Fe— 03 dz 0,548 
06 — Bi3 ds 0,009 Fe— 04 dg 0,544 
06 — Bi4 dá 0,050 Fe—O5 dg 0,545 
Fe-—O1 d5 0,572 Fe — 06 dto 0,555 


Tabela 5-18: Mínima Densidade Eletrônica (MED) calculada para a amostra (Bi, Nd, )Feo, 


Na Figura 5-16 vemos os mapas de densidade eletrônica para a amostra (Bi, Nd, ,) 
FeO,, com o aumento na substituição por neodímio uma nova distorção entre os íons 
acontece, podemos observar na figura 5-16(c) ligações covalentes entre os íons de oxigênio 
e bismuto. Nas figuras 5-16 (b) e (d) apesardas distorções é possível identificar as ligações 
de covalência entre os íons de ferro e oxigênio. 
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(a) (100) (b) (200) 


(c) (001) 


Figura 5-16: Mapas de contorno bidimensionais de densidade eletrônica da amostra (Bi, Na, )Feo,. 


Na figura 5-17 vemos variação da densidade eletrônica na linha de conexão dos 
íons adjacentes da amostra (Bi, ,Nd, )FeO, . As medidas de mínima densidade eletrônica 
são apresentadas na Tabela 5-19, se comparamos com as medidas da amostra anterior 
que possuía 10% de neodímio veremos que a única ligação que continua iônica é a d, a 
ligação dº é iônica com caráter covalente e as ligações d,e d, passaram a ser covalentes 
talvez pela presença dos Lone Pairs. As ligações entre o íon de ferro e os íons de oxigênio 
são covalentes, mas com menor intensidade se comparadas com a amostra (Bi, Nd,,) 
FeO 


3 
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“85 1.0 1.5 2.0 2.5 3.0 


d(Á) 
Figura 5-17: Variação da densidade eletrônica na linha de conexão dos íons adjacentes na amostra 
(BiyaNd, )JFeo.. 
Ligação Símbolo MED (CÂ) Ligação Símbolo MED (CÂ) 
06 — Bit [oh] 0,048 Fe— O2 d6 0,401 
06 — Bi2 do 0,135 Fe-— Os az 0,322 
06 — Bi3 ag 0,297 Fe — 04 dg 0,197 
06 — Bi4 d4 0,207 Fe— 05 dg 0,223 
(Ch ds 0,362 Fe— 06 d10 0,219 


Tabela 5-19:Mínima Densidade Eletrônica (MED) calculada para a amostra (Bi, jNd, ,)FeO.. 


Na Figura 5-18 podemos observar os mapas de densidade eletrônica para a amostra 
(BissNd,;)(Feç9900501)05, NO plano (001) (figura 5-18(c)) verificamos uma ampla nuvem 
eletrônica para o íon de oxigênio. Para os planos (200) e (002) (figura 5-18(b) e (d)) vemos 
ligações covalentes entre os íons de Fe*+ e O2. 
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(a) (100) (b) (200) 


Figura 5-18: Mapas de contorno bidimensionais de densidade eletrônica da amostra (Bi, Nd,,) 
(Fes 9009 01)05: 


0,99 


Na Tabela 5-20 observamos a mínima densidade eletrônica calculada para as 


» vemos que todas as ligações são 


ligações da amostra de (Bi, Nd, (Fes gsC0,0:)0 


0,99 0,01 
covalentes com interações maiores que em todas as outras amostras vistas neste trabalho. 
Na figura 5-19 observamos a variação da densidade eletrônica na linha de conexão dos 
O,. 


fons adjacentes da amostra (Bi, Nd, ,)(FesgsC0901) 
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2.0 


“8, 


[e 


1.0 


3.0 


Figura 5-19: Variação da densidade eletrônica na linha de conexão dos íons adjacentes na amostra 
(Bi, Nd, (Fes gs 0001)0s- 


Ligação Símbolo MED (CÁ) Ligação Símbolo MED (CÂ”) 
06-Bil d7 0,472 Fe —- O2 d6 1765) 
06- Biz d2 0,495 Fe-— 03 dz 1,135 
06- Bi3 ds 0,359 Fe-— 04 dg 0,728 
06 -— Bi4 d4 0,581 Fe-— 05 dg 0,734 
Fe-—O1 d5 1,154 Fe — 06 d10 0,754 


Tabela 5-20: Mínima Densidade Eletrônica (MED) calculada para a amostra (Bi, Nd, ,)(Fesgs00501)0s: 


Observando os valores de MED da tabela 5-21, e comparando-os com as distâncias 


de ligação da tabela 5-15, podemos notar que não há uma correlação evidente entre estas 


duas grandezas. Ou seja, uma distância menor do que aquelas apresentadas na tabela 


5-15 não necessariamente leva a uma ligação covalente, assim como distâncias maiores 


não necessariamente representam ligações iônicas. 
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Distâncias de ligação entre os íons do 


Ligação Composto (Bi, Nd,J(Fe, Co))O, (À) 
x=0,0e x=0,0e x=0,0:e x=0,1 e x=0,2e x=0,1 e 
y=0,0 y=0,01 y=0,02 y=0,0 y=0,0 y=0,01 
Bi—O 2,13 2,173 2,13 2,12 Ed 2,12 
Fe—O 2,04 2,04 2,04 2,04 2,04 2,04 


Tabela 5-21: Distâncias de ligação calculadas com base no raios iônicos de Shannon. 


q “5 
MED (€A ) para as Dopagens do Composto 


Ligação Símbolo (Bi, Nd (Fe, ,Co)JO, 
x=0,0ey= |x=00ey=| x=0,1ey x=0,2ey | x=0,1ey 

0,01 0,02 =(0)/0) = 0,0 = 0,01 
06 — Bi d1 0,087 0,155 0,005 0,048 0,472 
06 — Biz da 0,080 0,070 0,003 0,135 0,495 
06 — Bi3 dg 0,160 0,051 0,009 0,297 0,359 
06 — Bi4 d4 0,243 0,047 0,050 0,207 0,581 
Fe-—-O1 ds 0,493 0,914 0,572 0,362 1,154 
Fe — 02 de 0,475 0,914 0,581 0,401 iIRiiZo 
Fe-— O3 d7 0,470 0,890 0,548 0,322 1,135 
Fe— 04 dg 0,154 0,300 0,544 0,197 0,728 
Fe— 05 dg 0,154 0,310 0,545 0,223 0,734 
Fe — 06 d10 0,116 0,316 0,555 0,219 0,754 


Tabela 5-22: MED (eÃ*) para as Dopagens do Composto (Bi. Nd(Fe,,Co)O, 
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CONCLUSÕES 


Neste trabalho realizamos os estudos sobre as estruturas cristalinas do compostos 
(Bi, Nd (Fe, ,Co)O, com (x=0,0;0,1€e 0,2)e (y = 0,0; 0,01 e 0,02), utilizando difração de 
raios X e refinamento método de Rietveld. 

Comparando nosso trabalho com as propostas apresentadas na tabela 2- 3 para 
estruturas do composto (Bi, Nd,)FeO,, principalmente em relação as dopagens com 20% 
de neodímio, a qual as referências [23], [32]-[34], [36] indicaram as fases ortorrômbica e 
triclínica para representar a estrutura, verificamos que não houve mudança na estrutura do 
composto. Tendo permanecido a fase romboédrica (R3c) como a que melhor representa 
a estrutura com as dopagens para x = 0,1 e x = 0,2. Para o composto Bi(Fe,,Co))O, 
obtivemos o mesmo resultado que Luo e colaboradores [37], ou seja, o grupo espacial 
que melhor representa a estrutura do composto dopado com cobalto é o romboédrico 
(R30). Em relação as amostras com as dopagens simultâneas de neodímio e cobalto, Xue 
e colaboradores [38], mostraram que houve uma transição de fase para o composto (Bi, 
Nd (Fe, Co 


9970090904 de romboédrica para triclínica com o aumento do teor de neodímio, 


o mesmo não foi observado por nós, onde a amostra com x = 0,1 e y = 0,01 apresentou 
estrutura romboédrica (R3c) e a amostra com x = 0,2 e y = 0,02 não foi possível identificar, 
até o momento, qual estrutura melhor representa. 

Com relação a natureza de ligação entre os íons do composto (Bi, Nd)(Fe,,Co,) 
O, a análise foi realizada a partir dos cálculos de densidade eletrônica fornecidos por 
meio do programa Dysnomia. Ainda utilizamos os raios iônicos de Shannon [54], como 
critério para distinguir as ligações entre iônica e covalente. Pudemos observar que com a 
adição de cobalto as ligações entre os íons de ferro e oxigênio tendem a se tornar mais 
covalentes. O mesmo não acontece com as ligações entre os íons de bismuto e oxigênio. 
Para as amostras dopadas com neodímio, pudemos observar que com x = 0,1 houve a 
total supressão da ligação covalente (Lone Pairs) entre os íons de bismuto e oxigênio 
e um aumento significativo na covalência entre os íons de ferro e oxigênio, a simetria 
entre os íons, para x = 0,1, caracterizam uma estrutura cúbica. Com essa característica o 
material deixa de ser ferroelétrico e passar a ser paraelétrico, mas trabalhos futuros são 
necessários para comprovar se isso realmente acontece. 

As perspectivas para a continuação deste trabalho, será os estudos de refinamento 
por método de Rietveld das estruturas cristalinas dos compostos (Bi, Nd, ,)FeO, e (Bi, ;Nd, ») 
(Fe, ..Co 


098009 02)0;: e ainda um estudo de densidade eletrônica obtidos por MEM (experimental) 


e DFT (teórico) para um estudo mais preciso da estrutura cristalina do composto (Bi, Nd,) 
(Fe,,Co)O.. 
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